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ABSTRACT 

This  is  an  overall  study  of  the  whole  process  of  projectile  pene¬ 
tration  into  soil.  The  governing  differential  equations  for  the  soil 
target  medium  are  developed  using  the  equations  of  motion  and  the 
continuity  condition.  Constitutive  relationships  are  derived  from 
Mohr-Coulomb  theory  for  a  rigid  plastic  material  and  from  viscosity 
theory.  The  other  necessary  equations  are  generated  from  an 
experimental  investigation  of  the  penetration  mechanism. 
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SUMMARY 


Recognizing  that  aerospace  power  supplies  (both  reactor  and 
isotope  types)  can  be  subjected  to  earth  impact  as  a  result  of  circum¬ 
stances  that  may  occur  with  a  reasonable  degree  of  probability,  Sandia 
Corporation  early  became  concerned  with  the  study  of  earth  impact. 
More  recently,  with  the  increasing  interest  in  the  philosophy  of  intact 
reentry  and  disposal  by  penetration  into  the  earth's  surface,  Sandia' s 
investigations  were  amplified  to  include  earth  penetration. 

One  of  the  objectives  of  the  current  investigation  is  to  develop 
an  equation,  utilizing  dynamic  properties  of  earth  materials,  that 
will  predict,  with  a  reasonable  degree  of  accuracy,  the  depth  of  pen¬ 
etration  achieved  by  a  projectile  with  a  known  trajectory  impacting 
any  given  point  on  the  earth's  surface. 

This  report  describes  the  work  that  has  been  done  to  date  in 
the  development  of  an  analytical  expression  which  describes  the  earth 
penetration  process. 
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CHAPTER  I 


STATEMENT  OF  PROBLEM  AND  HISTORICAL  BACKGROUND 

I.  INTRODUCTION 

The  basic  problem  in  penetration  studies  is  to  determine 
whether  the  total  force  on  a  given  projectile  can  be  calculated  at  an 
instant  of  time  when  it  is  at  some  depth,  moving  at  a  given  velocity 
in  a  given  target  by  using  real  parameters  that  describe  the  mechan¬ 
ical  behavior  of  the  target  during  penetration. 

The  purpose  of  this  study  is  to  formulate  a  complete  set  of 
equations  that  might  realistically  describe  the  behavior  of  the  soil 
target  and  the  projectile  during  the  penetration  event  by  using  equations 
from  continuum  mechanics  and  soil  mechanics  and  by  studies  that  will 
delineate  the  mechanism  of  penetration  and  the  boundary  conditions  that 
those  equations  must  satisfy.  The  simplest  case  of  a  vertical  pene¬ 
tration  of  a  rigid  right  circular  cylinder  into  a  half  space  of  homo¬ 
geneous  isotropic  soil  was  used  in  this  study.  Some  consideration  will 
be  given  to  other  projectile  configurations.  The  projectile  is  con¬ 
sidered  to  be  a  rigid  body  moving  at  a  subsonic  velocity  with  respect 
to  the  target  medium  in  free  flight  with  its  cylindrical  axis  parallel 
to  the  direction  of  gravity. 
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Current  interest  in  soil  penetration  problems  stems  from  the 
following  types  of  questions: 

1.  How  thick  must  an  embankment  be  or  how  much ’soil  cover 
is  required  to  protect  forces  and  equipment  from  artillery 
fire? 

2.  What  forces  are  exerted  upon  a  projectile  as  it  impacts  into 
soil,  water,  and  soil  water  air  mixtures?  How  can  the 
projectile  be  designed  structurally  to  withstand  these  forces? 

3.  When  a  space  nuclear  power  supply  is  returned  from  orbit, 
what  will  be  the  depth  of  penetration?  How  can  the  pro¬ 
jectile  be  designed  structually  to  withstand  these  forces? 
so  as  to  keep  it  intact  and  achieve  the  required  depth  of 
penetration? 

4.  When  a  penetrator  impacts  into  an  unknown  soil  and  an 
acceleration-time  signature  is  transmitted,  is  it  possible 
to  determine  the  physical  properties  of  the  soil  and  make 
an  estimate  of  its  bearing  capacity  and  mobility  ? 

These  are  only  a  few  of  the  problems  that  can  be  posed,  but  these 
will  serve  to  illustrate  the  current  interest  in  the  phenomenon  of  pene¬ 
tration.  If  the  basic  problem  can  be  solved,  many  subsidiary 
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problems  can  be  resolved,  and  it  will  be  possible  to  predict  the 
terminal  trajectory  and  depth  of  penetration. 

Penetration  may  be  defined  as  the  forced  movement  of  an  object 
into  a  larger  object.  The  problem  may  be  classified  according  to  the 
amount  of  penetration  and  the  deformation  of  the  projectile. 

1.  The  projectile  may  hit  the  target  and  bounce  back,  con¬ 
serving  its  energy.  This  is  a  perfectly  elastic  impact  with 
no  deformation  of  target  or  projectile. 

2.  The  projectile  may  dent  the  target  or  form  a  small  crater, 
then  bounce  back,  and  the  projectile  may  be  deformed. 

3.  The  projectile  may  impact  and  continue  its  terminal  trajec¬ 
tory  until  all  its  energy  is  used  and  it  stops  or  comes  to 
terminal  velocity  as  in  a  fluid  target. 

The  problem  may  be  classified  according  to  the  forcing  function 
moving  the  projectile  at  the  time  of  impact. 

1.  The  forcing  function  G(t)  may  be  zero — free  flight  caused 
by  some  previous  input  of  energy. 

2.  The  G(t)  might  be  a  vibratory  load — in  the  case  of  driving 
a  pile  or  vibrating  a  footing  on  soil. 
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3,  The  G(t)  might  be  a  single  transient  load — such  as  an  object 
placed  on  a  target  and  hit  with  a  hammer  or  even  a  nuclear 
blast  shock  wave  passing  over  the  object. 

4.  The  G(t)  might  be  a  continuously  increasing  force  with  time — 
as  in  the  case  of  a  rocket  being  thrusted  at  the  time  of 
impact  or  in  the  case  of  a  man  forcing  a  spade  into  the 
ground. 

The  problem  may  be  classified  according  to  the  angle  between  the 
trajectory  and  the  direction  of  gravity,  the  angle  between  the  face  of 
the  target  and  the  direction  of  gravity  and  the  angle  between  the  trajec¬ 
tory  and  the  face  of  the  target.  The  mechanical  strength  of  the  target 
influences  the  mechanism  of  penetration  in  soils  because  the  strength  is 
stress  dependent.  Therefore,  the  configuration  of  the  target  and 
restraints  to  movement  at  the  boundary  of  the  target  influence  the  soil 
strength  and  the  mechanism  of  penetration.  For  example,  shooting 
a  projectile  into  the  side  of  a  bank  is  different  from  vertical  pene¬ 
tration  into  a  half  space  of  soil.  See  Figure  1. 

Thus  there  are  many  types  of  penetration,  and  in  some  instances 
the  mechanisms  of  penetration  are  similar.  However  the  earth  pene¬ 
tration  problem  is  unique  in  that  no  other  type  of  penetration  completely 
duplicates  it.  The  penetration  of  armor  plate  by  a  projectile  is  a 
shock  wave  phenomenon  and  a  different  problem  from  penetration  into 


Longitudinal 

elevation 


Symmetrical  Target  Restinj! 
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a  half  space.  Hypersonic  impact  of  a  meteor  into  metal  probably  is 
not  the  same  problem  as  impact  of  a  .  22  caliber  bullet  into  metal  even 
though  there  may  be  similarities  in  the  phenomena. 

The  mechanism  of  earth  penetration  is  extremely  complicated;  in 
the  past  it  has  received  little  attentiono  Even  for  the  simplest  case  of 
vertical  penetration  into  a  half  space,  the  velocity  regime  at  impact 
undoubtedly  plays  a  dominant  role  in  determining  the  mechanism  by 
which  the  soil  is  displaced  as  the  projectile  passes  from  air  into  soil. 
Consequently,  it  cannot  be  expected  that  a  single  comprehensive  theory 
can  be  readily  evolved  which  will  predict  the  penetration  process  in  its 
totality  for  all  materials,  impact  velocities,  and  collision  geometries.^ 

In  the  past,  most  investigators  have  concentrated  on  developing 
theoretical  relations.  Some  of  these  relationships  are  rational,  but 
most  of  them  are  empirical,  deduced  from  available  data  and  dimen¬ 
sional  analysis  and  presumed  to  apply  only  within  definite  limits  of 
initial  velocity.  Target  geometry  has  essentially  been  ignored.  A 
rational  description  of  the  phenomenon  of  penetration  can  not  be  ob¬ 
tained  if  only  the  depth  of  penetration  and  initial  velocity  are  measured 
and  the  target  perhaps  examined  after  the  event  has  occurred.  The 
divergent  results  obtained  by  different  investigators  from  experimental 

^Werner  Goldsmith,  Analytical  Versions  of  Penetration  Processes, 
NAVWEPS  Report  7812,  NOTS  TP  2811,  Febru^y  1962,  p.  1. 
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evidence  of  supposedly  virtually  identical  tests  are  due  to  lack  of 
knowledge  of  the  mechanism  of  penetration.  Also,  it  is  very  difficult 
to  build  identical  targets  of  soil  or  to  hit  the  same  target  area  in  full- 
scale  tests.  Recently  it  has  become  possible  to  measure  acceleration 
or  position  of  the  projectile  during  terminal  trajectory  in  soil,  but  the 
response  of  this  instrumentation  is  limited,  and  many  problems  are 
yet  to  be  solved  if  exact  records  are  to  be  obtained. 

In  this  study,  equations  were  formulated  in  an  attempt  to 
describe  realistically  the  behavior  of  the  target  during  the  pene¬ 
tration  event.  Two-dimensional  tests  were  performed  so  that  the 
penetration  mechanism  could  be  recorded  with  high-speed  photography. 
Later,  three-dimensional  tests  were  performed.  These  tests  led  to  a 
description  of  the  boundary  conditions  and  the  effects  of  target  geom¬ 
etry  and  target  constraints.  Equations  were  then  formulated  and 
assumptions  were  made  which  contribute  to  the  solutions  of  the 
equations. 

It  must  be  remembered  that  the  terminal  trajectory  penetration 
is  a  nonsteady-state  problem  with  one  force-free  boundary.  Even  in 
the  case  of  fluids,  few,  if  any,  nonsteady  problems  have  ever  been 


solved. 
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II.  HISTORICAL  BACKGROUND 
The  problem  of  earth  penetration  by  a  projectile  received 

2 

attention  as  early  as  1742  in  Robins,  New  Principles  of  Gunnery. 

The  great  mathematician  Euler  considered  the  problem  in  his 

Neue  Grundsatze  Per  Artillerie  in  the  year  1745.  Poncelet  studied 

2 

the  problem  in  1835;  Petry's  analysis  was  made  in  1910.  From 
1910  to  1943,  Peres,  Probert,  Morin  and  Didion,  H.  P.  Robertson, 
Nobile  de  Giorgi  and  others  reported  work  which  evaluated  the 

3 

constants  in  the  formulas  previously  developed.  All  were  concerned 
with  the  thickness  of  embankments  required  to  protect  facilities  and 
personnel  from  artillery  fire.  The  significant  findings  of  each  will 
be  summarized  in  the  paragraphs  which  follow. 

All  of  the  above  investigators  assumed  the  soil  target  to  be 
homogeneous  and  isotropic,  with  a  horizontal  surface  extending  to 
infinity  in  every  direction  and  infinitely  deep;  i.  e,  ,  a  homogeneous  soil 
half  space.  They  also  assumed  that  the  projectile  or  missile  impacts 

2 

H.  P.  Robertson,  Terminal  Ballistics,  National  Research 
Council,  2101  Constitution  Avenue,  Committee  on  Passive  Protection 
Against  Bombing,  Washington,  D.  C.  ,  January  1941. 

3^ 

Anonymous,  Terminal  Ballistics  and  Explosive  Effects,  National 
Research  Council,  Committee  on  Passive  Protection  Against  Bombing, 
Washington,  D.  C.  ,  Vol.  2,  Appendix  to  Final  Report,  1943. 
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normal  to  the  half  space  and  that  its  axis  is  tangent  to  the  trajectory. 
The  penetration  path  was  considered  to  be  vertical.  Tumbling  of  the 
projectile  was  not  considered.  The  missile's  shape  was  assumed  to 
remain  undistorted. 

Newton's  equation  of  motion  is  the  starting  point  for  all  these 
theories, 

M  g  -  F  =  M  a  (I-l) 

P 

which  must  hold  when  the  projectile  is  at  any  depth  and  moving  at  any 
velocity,  where 

M  =  mass  of  projectile  (M  )  plus  the  soil  moving  with  the 

P 

projectile  (M^) 

g  =  acceleration  due  to  gravity 
F  =  soil  resistance  to  movement 
a^  =  acceleration  of  projectile. 

The  free  body  for  this  equation  is  shown  in  Figure  2. 

Gravity  forces  on  the  projectile  are  small  compared  to  the  soil 
resistance  force;  thus.  Equation  I-l  was  rewritten 

d  V 

-F  =  M  V 

p  p  dz 


(1-2) 
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Figure  2.  Dynamic  Free  Body  of  Projectile  After  It  Has  Entered  the  Soil 
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where  is  the  velocity  of  the  projectile  and  z  is  the  distance  meas¬ 
ured  from  the  surface  downward  along  the  trajectory.  If  no  soil  is 
considered  to  be  moving  with  the  projectile,  M  becomes  or  the  mass 
of  the  projectile  only. 

Euler,  Robins  and  Peres  assumed  that  the  soil  resistance  or  drag 
on  the  projectile  was  independent  of  velocity  and  depth  at  any  instant  of 
time  and  proportional  to  the  maximum  Cross-sectional  area  of  the 
projectile.  Therefore,  F  =  C  A  where  C  =  a  constant  and  A  =  maximum 
area  of  the  projectile.  Equation  1-2  can  then  be  rewritten 

-  C  A  dz  =  M  V  d  V  (1-3) 

P  P  P 

when  z  =  0,  'the  velocity  at  impact. 

When  z  =  d,  V  =  0  where  d  is  the  depth  of  penetration.  Integrating 
Equation  1-3 


one  obtains 
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-  C  A  d  =  M 
P 


V 


2  .0 


1 


or 


C  A  d 


M 

P 


(1-4) 


Therefore,  in  the  Euler  equation  the  depth  of  penetration  is 


d  = 


W 

P 

2ACg 


(1-5) 


where  is  the  weight  of  the  projectile 

The  Euler  equation  assumes  a  constant  deceleration  (a  linear 

relationship  between  velocity  and  time)  and  a  parabolic  relationship 

between  di^lacement  of  the  projectile  and  time.  Since  the  soil 

resistance  force  is  assumed  to  be  independent  of  velocity,  it  must  be 

equal  to  the  static  resistance  to  movement.  But  this  static  resistive 

force  is  now  known  to  be  a  nonlinear  function  of  depth  and  area  of  the 

projectile.  If  d  is  measured  in  feet,  W  in  pounds,  A  in  square 

P 

2  2 
inches,  g  in  ft/sec  and  in  ft/sec,  Euler's  C  is  1900  Ib/in.  . 

Using  data  from  French  and  American  artillery  units  in  1932, 

2 

gave  C  as  2100.  0  Ib/in.  . 


Peres 
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Poncelet  assumed  that  the  soil  resistance  force  P  =  f.  (z)  f„(V  ) 

X  Z 

where  f  (V  )  is  a  function  only  of  the  velocity  of  the  projectile  at  any 
2  P 

instant  of  time  and  f  j^(z)  is  a  function  only  of  depth  at  the  same  time. 
Reasoning  that  the  fluid  drag  on  a  body  is  proportional  to  the  square 
of  the  velocity,  he  assumed 


f  (V  )  =  a  +  b  V 

2  p  p 

f^(z)  =  3  A 


(1-6) 


where  a,  b,  and  j  are  constants.  Substituting  these  relationships 
into  Equation  1-2,  one  obtains 


-jA(a  +  bV  )dz  =  M  V  dV 
'*  p  P  P  P 


V 


-  3  A  dz  =  M 


a  +  b  V 


d  V 

2  p 


P 


Integrating  as  before,  using  the  boundary  conditions  z  =  0  when 

V  =  V  at  impact  and  z  =  d  when  V  =0  and  when  the  total  depth  of 
pi  P 

penetration  is  reached  and  the  projectile  is  at  rest, 

2 

/  h 

j  A  d  =  Mp  ^ 


or 


W 


d  = 


2  Ajb 


jtn 


(^) 


(1-7) 


Probert,  Morin  and  Didion  tried  to  evaluate  the  a  and  b  param¬ 
eters  of  Poncelet's  equation.  The  j  was  ignored  or  considered  to  be 
equal  to  one.  H.  P.  Robertson  reviewed  their  data  and  decided  that 

V  ^ 

f(V  )  =  a  +  P-  (1-8) 

P  2  g 


where  k  is  a  new  constant  and  7  is  the  unit  weight  of  the  soil.  Thus 
Poncelet’s 


b 


7k 

2g 


and  Robertson's  formula  is 


A7k  \  2g  1  / 


(1-9) 


W 


In  1910,  Petry  called  =  o’,  the  frontal  loading  of  the  projectile 
in  pounds  per  square  inch,  and  changed  the  formula  somewhat  to  obtain 

2 

^  (I-IO) 
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If  a  is  expressed  in  pounds  per  square  inch,  d  in  feet  and  in 

2 

ft/sec.  ,  thert  the  k  is  in  ft.  in.  /pound.  Values  given  for  k  were 


5.  30  for  a  "sandy  soil" 

6.  95  for  a  "soil  containing  vegetation" 
10.  6  for  a  "clay  soil" 


0.  775  for  limestone 
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The  Petry  formula  is  simple.  A  nomogram  can  be  prepared 
relating 

a.  frontal  loading  (weight  of  projectile  /frontal  area) 

b.  impact  velocity 

c.  depth  of  penetration 

d.  constant  of  proportionality 

The  procedure  was  to  perform  a  test  measuring  the  first  three 
parameters,  then  to  calculate  "k."  The  k  is  a  variable  "constant"  which 
depends  on  many  parameters.  Therefore,  to  make  a  prediction  of  the 
depth  of  penetration,  it  is  almost  necessary  to  have  performed  a 
similar  test  before,  using  the  same  projectile,  soil  and  impact  velocity. 
In  the  absence  of  a  formula  based  on  theory  and  utilizing  the  proper 
soil  parameters,  prediction  will  still  have  to  be  made  this  way. 

This  is  the  empirical  approach. 

In  summary,  for  all  of  these  penetration  formulas  it  is  assumed 
that  the  soil  resistance  is  independent  of  depth  and  directly  proportional 
to  the  area  of  the  projectile.  In  Euler's  formula,  the  soil  resistance 
is  independent  of  velocity;  and  in  all  the  others,  the  soil  resistance 
force  is  a  linear  function  of  the  square  of  the  velocity.  The  Euler 
formula  contains  one  parameter  which  is  independent  of  velocity; 
therefore,  it  must  be  equal  to  the  static  resistance.  But  it  is  now 
known  that  the  static  resistance  of  soil  is  a  nonlinear  function  of  area 
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and  depth.  It  is  even  doubtful  that  there  is  any  simple  relationship 
between  static  resistance  and  dynamic  resistance,  even  if  inertial 
resistance  of  the  soil  is  taken  into  account.  In  the  "static"  case,  the 
soil  is  essentially  at  rest  or  moving  with  a  very  small  velocity.  In 
the  dynamic  case,  the  soil  particles^  have  considerable  velocity  which 
can  affect  the  mechanism  of  penetration. 

In  the  Poncelet  formula  there  are  two  constants.  The  "a"  is 
considered  to  be  the  static  resistance  independent  of  velocity  and 
depth,  and  the  "b"  a  multiplier  for  the  velocity  square  term.  The 
same  objections  would  hold  true  for  the  Poncelet  "a"  as  were  given 
above  for  the  Euler  constant.  The  drag  on  a  body  moving  through  a 
fluid  is  a  linear  function  of  the  velocity  squared;  but  this  applies  to  the 
steady- state  case  only  where  the  velocity  is  constant,  where  the  body 
is  submerged  in  a  full  space,  and  where  target  body  forces  are  un¬ 
important.  These  are  certainly  not  the  conditions  in  the  penetration 
event  where  the  projectile  is  slowing  down  in  a  half  space  and  the 
target  body  forces  due  to  gravity  probably  make  one  of  the  most 
important  contributions  to  the  resistance,  since  the  soil's  strength  is 
stress  dependent. 

In  all  of  these  earlier  attempts,  a  solution  for  the  soil  force 
resisting  movement  has  been  assumed  for  use  in  the  one  differential 
equation  of  motion  of  the  projectile.  This  equation  is  integrated 
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between  the  limits  of  initial  velocity  at  zero  depth  and  zero  velocity 
at  final  depth.  By  performing  tests  where  the  initial  velocity  and 
final  depth  could  be  measured  when  the  projectile  properties  were 
known,  it  was  possible  to  compute  the  "constants”  which  satisfied 
the  equation.  If  the  assumed  solution  is  correct,  it  should  hold  for 
all  soils,  projectiles,  and  initial  velocities. 

This  was  the  only  possible  approach  to  the  problem  in  the  past 
since  (1)  velocities,  displacements  and  accelerations  could  not  be 
measured  during  the  penetration  event  which  would  immediately  allow 
the  checking  of  the  assumed  force  solution;  (2)  there  was  no  knowledge 
of  the  mechanism  of  penetration — i.  e.  ,  how  the  soil  was  displaced 
from  in  front  of  the  projectile  and  how  this  mechanism  was  affected  by 
the  physical  boundary  restraints  on  the  target  and  target  geometry; 
and  (3)  there  were  no  constitutive  equations  to  describe  the  relationship 
between  forces  imposed  by  the  projectile  and  the  amount  of  deformation 
and/or  the  rate  of  deformation  of  the  soil  medium. 


CHAPTER  II 


THE  ANALYTICAL  APPROACH 

1.  INTRODUCTION 

Numerous  attempts  have  been  made  to  solve  the  basic  earth  pene¬ 
tration  problem  by  using  the  one  equation  of  motion  of  the  projectile, 
Equation  I-l.  There  are  three  unknowns  in  this  equation:  F  the  force 

on  the  projectile,  V  the  velocity  of  the  projectile  and  M  the  mass  of 

p  c 

the  soil  that  moves  with  the  projectile.  Thus,  two  additional  equations 
are  needed. 

The  static  (very  slow)  resistance  to  penetration  of  a  rigid  body 
such  as  a  pile  or  footing  into  a  half  space  has  never  been  rigorously 
solved.  However  good  approximations  are  available  for  shallow  pene¬ 
trations  and  were  developed  primarily  by  Prandtl  and  extended  by 
4 

Terzaghi  and  others.  In  this  slow  steady- state  solution,  the  force  P 
required  to  move  a  circular  plate  slowly  downward  is  given  by 

P=;rR^[1.3CN  +  7d  N  +0.  SyRN]  (II-l) 

c  q  y 

4 

Karl  Terzaghi,  Theoretical  Soil  Mechanics,  (New  York: 

John  Wiley  &  Sons,  1943),  Chap.  VIII. 
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where  N  ,  N  and  Ny  are  nonlinear  functions  of  the  angle  of  internal 

5 

friction  only 

C  is  the  cohesive  strength  of  the  soil 
d  is  the  depth  of  soil  measured  from  the  soil  surface 
to  the  bottom  of  the  loaded  plate 
y  is  the  effective  unit  weight  of  the  soil 
R  is  the  radius  of  the  loaded  circular  plate  and 
cp>  is  the  angle  of  internal  friction  of  the  soil. 

The  bearing  capacity  factors  N^,  and  are  shown  in  Figure  3. 

In  light  of  the  complexity  of  even  the  static  case,  it  seems  hope¬ 
less  to  assume  the  soil  resistance  function  F  for  even  one  set  of  con¬ 
ditions  in  the  dynamic  case.  The  solution  of  the  dynamic  penetration 
problem  requires  the  same  general  type  of  analysis  that  Prandtl  used 
for  the  static  case.  The  problem,  however,  is  much  more  complicated 
because  inertial  and  viscous  effects  must  be  considered.  One  must 
recognize  that  the  dynamic  penetration  problem  involves  the  nonsteady 
case.  The  projectile  is  slowing  down  in  proportion  to  the  force  exerted 
on  the  projectile. 

If  high  response  accelerometers  can  be  devised,  it  may  be  possi¬ 
ble  to  obtain  the  trends  in  the  function  F  for  various  weights  of  projectiles, 

5 

G.  G.  Meyerhof,  "The  Ultimate  Bearing  Capacity  of  Foundations," 
Geotechnique,  II  (1951),  p.  301. 
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areas,  nose  shapes,  entrance  velocities  into  various  targets  with  dif¬ 
ferent  geometries  and  mechanical  properties.  These  accelerometers 
are  just  now  being  developed,  but  still  the  amount  of  testing  required 
for  an  empirical  solution  to  the  basic  problem  will  be  monumental. 
Therefore,  it  seemed  reasonable  to  use  available  continuum  equations 
for  the  target  medium,  and  constitutive  equations  for  soil  media  to  gene¬ 
rate  equations,  the  solution  of  which  may  be  helpful  in  the  reduction 
of  the  monumental  testing  program  otherwise  required. 

II.  GENERAL  CONTINUUM  EQUATIONS 
Equations  of  Motion  of  the  Soil  Target  Medium 

First,  consider  Newton's  second  law  written  for  the  soil  particles 

for  the  three  spatial  orthogonal  directions.  The  simple  and  elegant 

equation  of  motion  for  any  continuous  medium  was  first  discovered  by 
0 

Cauchy.  It  is  given  by 

dV  —  = 

p-^  =  pB  +  div  T  (II-2) 

where  p  =  the  soil  density  function 

V  =  the  velocity  vector 
B  =  the  body  force  vector 

^A.  -L.  Cauchy;  Ex.  de  Math.  3  (1828),  Oeuvres  (2)  8, 
pp.  195-226. 
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div  =  the  divergence  operator 
T  =  the  stress  tensor 


This  equation  is  valid  independent  of  the  mechanical  properties, 
of  the  material  and  the  form  the  stress  tensor  may  take. 


Cylindrical  coordinates  have  been  adopted  for  the  inertial  space, 
and  the  Eulerian  or  spatial  form  employed  is  shown  in  Figure  4.  The 
body  force  vector  is  equal  to  the  acceleration  due  to  gravity  acting 
downward  in  the  positive  z  direction  or  parallel  to  the  trajectory  of  the 
missile.  Thus  Cauchy's  equation  in  the  three  cardinal  directions 
becomes  the  commonly  known  equations  of  motion  in  cylindrical 
coordinates 


(II- 3) 


(11-4) 


(II- 5) 


where  compressive  stresses  are  taken  as  positive.  The  equations  are 
valid  for  any  point  in  the  target  space  where 


"^J.  S.  Sokolnikoff, 
McGraw-Hill  Co.  ,  1956), 


Mathematical  Theory  of  Elasticity  (New  York: 
p.  184. 


Gravity  in  z  direction 
t^  =  time  of  surface  impact 
t  =  any  time  after  t 

d  =  J  Vp  dt  “  velocity  of  projectile 

^o 

z  =  Z  -  d 


Figure  4.  Cylindrical  Coordinate  System  for  the  Soil  Space 
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o  - 

r 


a  = 
■  z 


V  = 
r 

V  = 

V  = 
z 

P  = 


g  = 


normal  stress  in  the  r  direction 
normal  stress  in  the  ij*  direction 
normal  stress  in  the  z  direction 
velocity  component  in  the  r  direction 
velocity  component  in  the  ijj  direction 
velocity  component  in  the  z  direction 
density  of  the  medium  which  is  function  of  space 
and  time 

acceleration  due  to  gravity 


It  is  specifically  assumed  that  the  stress  tensor  is  symmetric; 
this  is  true  when  no  extraneous  couples  are  present  as  in  the  case  of  a 
polarized  medium.  Thus, 


T  =  T  =  shear  stress  in  the  r-z  plane 
rz  zr 


r  =  T  =  shear  stress  in  the  r-  iJ;  plane 

r  ,  =  T  =  shear  stress  in  the  z-4>  plane 

4JZ  zif 


The  case  considered  here  is  axisymmetric,  thus  t 

O 

and  =  0,  and  any  derivative  with  respect  to  lA  is  zero. 
Thus,  Equation  II-3  reduces  to 


=  0, 


(II-6) 


S.  Timoshenko  and  J.  N.  Goodier,  Theory _of  Elasticity, 
(New  York:  McGraw-Hill  Book  Co.  ,  1951),  p.  343. 
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Equation  II-4  reduces  to  the  identity  0=0. 
Equation  II-5  reduces  to 


(n-7) 


The  acceleration  terms  in  the  two  remaining  equations  of  motion 


.  9 


can  be  expanded  by  the  chain  rule  to  obtain 


dV  8  V 
r  r 


8  V 


8  V 


+  V 


+  V 


dt  8t  r  8r  z  8z 


(II- 8) 


d  V  8  V 
_ z  _ _ z 

dt  8t 


8  V 


8  V 


+  V 


+  V 


r  8r  z  8z 


(n-9) 


since  V  and  V  are  both  functions  of  (r,  z,  t). 
r  z 


Therefore,  the  equations  of  motion  can  now  be  rewritten  as 


da  dr 
r  rz 


1 

r  r  lu  '  dt 


8  T  da  . 

rz  z  1  ,  „ 

- T  +  pg 

8r  8z  r  rz 


=  (■ 


fdv  8  V 

8V  ' 

T  V 

1  +  V 

+  V  -  ^ 

8t  r  8r 

z  8z 

'dV  8  V 

8  V 

,  ^  +  V  ^ 

+  V  " 

8t  r  8r  z  8z 


Ip  (II- 10) 


p  (II- 11) 


Salamon  Eskinazi,  Principles  of  Fluid  Mechanics,  (Boston: 
Allyn  and  Bacon,  Inc.  ,  1962),  p.  143. 
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Continuity  Equation 

The  soil  possesses  a  density  function  p  =  p(r,  z,  t)  which  is  given 
by  the  formula 

M  =  I  p  d  V  (11-12) 

•'v 

where  M  =  the  mass  of  the  soil  particles  enclosed  within  a 

certain  region  v  with  the  physical  dimensions  of 
mass  per  unit  volume 

V  =  the  volume  enclosing  these  particles  and  is  not 
necessarily  a  constant 


p  =  mass  density  of  the  soil — not  necessarily  a  constant 
during  the  event 

The  physical  significance  of  the  concept  of  mass  is  given  by  the 

10 

principle  of  conservation  of  mass:  the  mass  of  the  soil  particles  in  a 
material  volume  v  does  not  change  as  v  moves  with  the  soil  particles 
being  displaced  by  the  projectile.  Thus  the  principle  of  conservation 
of  mass  is  otherwise  stated  by  the  expression 


(11-13) 


^*^James  Serrin,  "Mathematical  Principles  of  Classical  Fluid 
Mechanics,"  Handbuch  der  Physik,  (Berlin:  Springer  Verlag,  1959), 
Vin,  pp.  125-263. 
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Since  the  transport  theorem  is  valid ^ 

V  V 

where  G  is  any  scalar  or  vector  function 

V  is  the  vector  displacement  function 
div  is  the  divergence  operator 

Equation  11-13  becomes 


/( 


dt 


+  p  div  V  )  dr  =  0 


Since  v  is  arbitrary,  this  implies 


dt 


+  p  div  V  =  0 


in  Euler ian  or  spatial  form  and 


1^  +  div  (p  V)  =  0 


(n-15) 


(n-16) 


(n-17) 


in  Lagrangian  or  material  form. 

Using  the  Eulerian  form  and  cylindrical  coordinates,  the  con¬ 
tinuity  equation  becomes 


4e+  p( 

dt 


8  V  V  8  V 

8r  r  8z 


)=0 


(n-18) 


and  this  equation  is  valid  for  any  continuous  medium. 
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If  p  =  a  constant,  the  equation  becomes 


a  V  V  a  V 

+  —  +  — ^ 
ar  r  az 


=  0 


(11-19) 


The  derivation  given  above  is  substantially  due  to  Euler. 

The  Equation  of  Motion  of  the  Projectile  After  the  Projectile 
Has  Entered  the  Soil 

d  V 

(M  +M)g-F  =  (M  +M)  — (Equation  I-l)  is  the  equation  of 
pc  p  c  dt 

motion  of  the  projectile  as  it  moves  through  the  soil  where  F  is  the  sum 

of  all  forces  exerted  by  the  soil  upward  on  the  projectile.  This  force  F 

must  also  be  the  sum  of  all  forces  exerted  by  the  projectile  downward 

on  the  soil  except  the  weight.  The  term  M  is  the  mass  of  the  projectile, 

P 

and  is.mass  of  soil  that  might  be  carried  along  with  the  projectile.  But 


F  =  I  pda 


(11-20) 


where  is  the  closed  surface  of  the  projectile  and  any  soil  carried 
along  with  it,  da  is  the  differential  area  of  the  closed  surface,  and  p  is 
the  stress  vector  which  depends,  at  any  given  time,  on  the  position  and 
orientation  of  the  surface  element  da  making  up^.  If  n  denotes  the 
outward  directed  normal  to  j^,  then  p  =  p  (r,  z,  t;  n).  Thus,  if  the  stress 
field  is  known  at  any  instant,  then  p  can  be  calculated  if  ^  is  known. 


L.  Euler,  Principes  General  due  Movement  des  Fluids,  Hist. 
Acad.  Berlin  1755  (opera  omnia  II  12),  pp,  54-92. 
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Energy  Equation  or  Equation  of  State 

Most  of  the  work  done  in  deforming  the  soil  medium  is  trans¬ 
formed  into  heat,  and  it  is  recognized  that  there  will  be  an  increase  in 
temperature  of  the  soil  medium  during  the  penetration  process.  It  is 
assumed  that  thermal  stresses  are  small  in  comparison  with  stresses 
producing  plastic  flow,  and  that  temperature  changes  do  not  affect  the 
density  or  mechanical  properties  of  the  soil.  Interaction  between 
mechanical  and  thermal  processes  will  not  be  considered  in  this  study. 

Thus,  the  energy  equation  or  an  equation  of  state  will  not  be  used  as 

12 

an  additional  equation. 

Summary  of  Equations  Available  From  Continuum  Mechanics  and 
Requirements  for  a  Solution 

From  continuum  mechanics,  one  obtains  the  following  equations 
1.  Equation  of  motion  of  the  projectile 

dV 

(M  +M)g-F  =  (M  +M) 
p  c  °  p  c  dt 

(Equation  I-l) 

12 

William  Prager,  Introduction  to  Mechanics  of  Continua, 

(New  York:  Ginn  and  Company,  1961),  pp.  87,88. 
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2.  The  vertical  resisting  force  on  the  projectile  is  the  sum  of 
the  forces  exerted  by  the  projectile  on  the  soil 


(Equation  11-20) 


where  p  =  p  (r^  z,t;  n)  for  the  axisymmetric  case. 

3.  The  two  equations  of  motion  of  the  target  medium,  which 
are  independent  of  target  mechanical  properties 


da  dr  (a  -  a  , )  /O  V 
_ r  rz  r  ’  4*  _  f  ^ 

r  \ 


3  V 


9  V 


9r  9z 
(Equation  II- 10) 


9t 


+  V 


r  9r 


+  V 


z  9z 


dr  da  r  / 

,  rz  z  rz  ,  I 


(Equation  II- 11) 


9  V 
_ 2 

9t 


9  V 


9V 


+  V 


r  9r 


+  V 


z  9z 


4.  The  continuity  equation,  which  also  is  independent  of  target 
mechanical  properties 


,  /9V  V 

p  — 
dt  '  3r 


r  r 
+  —  + 


9  V 

2 

9z 


0 


(Equation  11-18) 


The  independent  variables  are  r,  z,  t.  The  dependent  variables  are 

M,F,  V,(T,T  ,0’,cr|,V,V  and  p.  Initially,  there  are 
c'  p'  r^  rz'  z'  vIj'  r'  z  ^ 

available  five  independent  equations  with  ten  unknowns.  It  is  also 
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required  that  the  surface  be  described  as  well  as  the  initial  con¬ 
ditions  (depth  of  penetration  zero)  of  velocity  and  stress.  The  boundary 
conditions  on  the  stress  and  velocity  fields  must  also  be  available  if  F 
is  to  be  calculated  or  the  basic  problem  solved. 

Thus,  five  additional  equations  are  required.  These  equations 
must  be  obtained  from  constitutive  relationships  and  assumed  relation¬ 
ships  between  the  unknown  parameters.  The  constitutive  equations 
relate  stresses  to  strain  and  strain  rate.  By  taking  advantage  of  the 
body  of  knowledge  relating  stress  to  yield  strain  in  soil  mechanics  and 
using  theory  from  rheology,  constitutive  equations  can  be  postulated. 
The  additional  required  equations  and  boundary  conditions  must  be 
developed  and  justified  in  the  mechanism  of  penetration  studies. 

III.  CONSTITUTIVE  RELATIONSHIPS  FOR  SOIL 
State  of  Stress  on  an  Oblique  Line  in  the  r-z  Plane 

By  equating  forces  on  an  element  of  soil  using  cylindrical 
coordinates  as  shown  in  Figure  5,  it  is  possible  to  determine  the 
normal  stress  cr^  and  the  shear  stress  on  the  line  inclined  to  the 
z  axis  at  angle  ^  in  terms  of  t  ,  a  ,  a  and  the  angle  /3.  As  shown 
by  calculations  on  page  33 ,  the  stresses  on  an  oblique  line  in  the 
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Consider  Case  I  of  Figure  5. 

Area  of  top  surface  of  element  =  r  dr  dijj 


Area  of  vertical  surface  of  element  =  rdij/dz 


Area  of  side  of  element  = 


dr  dz 
2 


Area  of  oblique  surface  = 


r  dr  dip  _  r  dll'  dz 
sir^TT-0j  cos|:7-/3^ 


Equating  Forces  in  r  Direction 

a  r  d4<  dz  +  T  r  dr  d4<  +  2  <t 


/dr  dilj  \  . 

2  / 


difj 


■  "^rz  ^  tan  ^  =  0 


since  the  term  is  of  higher  order 
Equating  Forces  in  z  Direction 


or^  r  di  dip  +  r  dip  dz  + 


r  dip  dz  cos  {n  -  /f) 
cosp:  -(3) 


r  dll'  dr  .  /  ,,\  _ 

-  a  ■■  a  smiff  -  fl)  =  0 
s  sin^-/^  '  '  ' 


CT  -  T  cot  B  -  T  cot  /3  -  cr  =0 
z  rz  s  s 


Solving  Simultaneously 


XT  +  O  /  XT  -  <T  \ 


s  2 


cos  2/5  -  T  sin  2/3 
rz 


XT  -  a 

T  =  ■  sin  2  /3  -  T  cos  2  B 

s  2  rz 


Figure  5  (cont).  Calculations  for  the  State  of  Stress  on  an  Oblique 

Plane  in  the  Axisymmetric  Cylindrical  Coordinates 
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r-z  plane  do  not  depend  on  the  azimuthal  normal  stress  for  the 
axisymmetric  case.  Thus  the  equations 

CT  +  (T  a  -  a 

— -+  ^2  ^  cos  2/6  -  sin  2;3  (11-21) 

a  -  a 

Sin  2  /3  -  T  cos  2iS  (11-22) 

S  ^  X'Z 

are  exactly  the  same  for  ct  and  t  in  the  two  dimensional  plane  strain 

s  s 

13 

state  of  stress.  Thus  the  two-dimensional  Mohr  circle  can  be  used  to 
represent  the  state  of  stress  if  it  is  kept  in  mind  that  cr ,  can  exist  but 
does  not  influence  the  normal  and  shear  stresses  on  the  oblique  line  in 
the  r-z  plane.  The  Mohr  circle  representation  of  this  state  of  stress 
is  shown  in  Figure  6.  These  equations  are  independent  of  the  material's 
mechanical  properties. 

Mohr-Coulomb  Yield  Condition 

The  Mohr-Coulomb  assumption  which  is  used  to  explain  the 
behavior  of  soil  at  failure  has  been  repeatedly  verified  for  so-called 
"static"  equilibrium  states  of  stress.  Essentially  this  is  a  stress- 
dependent  "law"  giving  the  combination  of  shear  and  normal  stress  for 
which  the  soil  yields  plastically.  The  formula  is  a  linear  representation 

13 

A.  R.  Jumikis,  Mechanics  of  Soils,  (Princeton,  N.  J.  : 

D.  Van  Nostrand  Co. ,  1964),  p.  14. 


Figure  6.  Mohr  Circle  Representation  of  State  of  Stress  in  r-z  Plane 
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of  the  failure  conditions  and  is  normally  satisfactory  for  practical  work 
where  failure  by  shear  or  rupture  is  involved. 

T  =  C  +  or  tan^  (11-23) 

s  s 

where  t  =  shear  stress  on  the  surface  of  failure  (sometimes 

s 

called  the  slip,  yield  or  failure  surface)  at  the  time 
of  failure 

C  =  cohesive  strength  which  is  dependent  on  the  past 
stress  (history  and  moisture  content  of  the  soil 

cr  =  normal  effective  stress  on  the  surface  of  failure  at 
s 

time  of  failure 


where  =  total  normal  stress  on  the  failure  or  yield  surface 

=  force  pressure  in  the  fluid  filling  the  voids  within  the 
soil  at  the  point  in  question 

and 

0  =  angle  of  internal  friction  (nearly  constant  for  a  given 
soil  with  a  given  density) 

In  Figure  7a  some  typical  plots  of  the  shearing  resistance  versus 
shearing  deformation  are  shown.  Henceforth  all  stresses  will  be  con¬ 
sidered  as  effective  stresses. 


Triaxial 

Test 


LX 

Dense  soil 


Direct  Shear  U 
Test 


Simple  Shear 
Test 


/^ultimate 


Loose  soil 

1 

1 

Normal  stress  on  failure  surface 
kept  constant . 

1 

Deformation 

c  a  for  simple  shear  test 

for  triaxial  test 
for  direct  shear  test 

a.  Typical  Stress -Deformation  Relationships  for  Soils 


Dense  soil 


Soil  at  critical  density 


1 

or  any  saturated  soil 

1 

1 

Loose  soil 

with  no  drainage. 

c  Deformation 

b.  Typical  Density-Deformation  Relationships  for  Soils 


Figure  7.  Typical  Density-  Stress- Deformation  Relationships  for  Soils 
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If  the  soil  is  dense,  its  peak  strength  is  greater  than  its  ultimate 
14 

strength.  Since  large  deformation  probably  occurs  during  the  pene¬ 
tration  event,  it  will  be  the  ultimate  strength  that  is  the  important 
parameter  and  should  be  used  in  determining  the  angle  of  internal 
friction.  By  performing  several  tests  and  then  plotting  the  shearing 
strength  t  ultimate  versus  the  effective  normal  stress  on  the  plane 
of  failure  at  the  time  of  failure,  the  Mohr-Coulomb  failure  envelope  is 
obtained  and  C  and  ^  ultimate  can  be  determined.  A  typical  plot  is 
shown  in  Figure  8.  Typical  plots  of  density  versus  deformation  are 
shown  in  Figure  7b  and  indicate  that  volume  changes  normally  occur 
during  at  least  the  early  phase  of  the  shearing  process.  If  the  stress 
given  by  Mohr-Coulomb  failure  criterion  is  taken  as  the  yield  stress 
of  an  ideally  plastic  medium,  the  nonhorizontal  part  of  stress- 
deformation  curve  of  Figure  7a  is  ignored.  The  deformation 

required  to  reach  ultimate  strength  is  normally  small  so  the  assump- 

15 

tion  above  seems  justified.  The  deformation  which  develops  in  the 
nonhorizontal  part  of  the  shear  resistance  deformation  curve  is  also 
usually  nonrecoverable  when  the  load  is  removed. 

^ ^Donald  Taylor,  Fundamentals  of  Soil  Mechanics,  (New  York: 
John  Wiley  &  Sons,  1948),  pp.  329-338. 

C.  Drucker  and  W.  Prager,  "Soil  Mechanics  and  Plastic 
Analysis  or  Limit  Design,"  Q.  Appl.  Math.  ,  X  (1952),  p.  157. 


Figure  8.  Yield  Line  Plotted  in  Mohr  Circle  Frame 
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The  Mohr-Coulomb  yield  condition  Equation  11-23  can  be  com¬ 


bined  with  Equations  11-21  and  11-22  if  the  stresses  ^  and  a  are  taken 

s  s 

as  the  stresses  on  the  yield  or  failure  surface;  thus  the  angle  I3 


designates  the  direction  of  the  yield  surface. 


From  Figure  9  it  can  be  seen  that  the  Mohr  stress  circle  repre¬ 


senting  the  state  of  stress  at  a  point  must  always  touch  the  line  of 


limiting  equilibrium  given  by  Equation  11-23  if  yielding  is  to  occur. 


If  the  stresses  on  the  bottom  side  of  the  element  are  always  considered 


positive,  i.  e.  ,  the  yielding  is  always  in  the  positive  r  direction,  then 


the  yield  line  must  pass  through  the  lower  point  of  tangency.  From 


Figure  9  it  can  be  seen  that,  at  yield. 


<7  =  CT  +  T  tan  6  +  r  cos  a 

z  s  s  max 

y 

(11-24) 

CT  =  CT  +  T  tand  -  T  cos  a 

r  s  s  max 

y 

(n-25) 

T  -  T  sin  a. 

rz  max 

y 

(11-26) 

But 

T 

s 

= - 1  ^ 

max  cos  <p 

(11-27) 

from  Equation  11-22 

T  =  C  +  CT  tand 
s  s 
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and 

c  =  2/3  +  ^-  -| 


(11-28) 


Then 


a  =0"  +  (C  +  CT  tan^) 
z  s  s 

y 


2  cos  j3  sin(^  -h  (j)) 
cos  0 


CT  = 

r 


y 


CT 

s 


(C  +  CT  tan  <j>) 
s 


2  sin  j8  cos(j3  + 
cos  ^ 


(11-29) 

(n-30) 


and 


T  =  -  (C  + 
rz 

y 


CT  tan  6) 
s 


cos  (2  /3  +  (^) 
cos  (^1 


(11-31) 


Thus  CT  ,  CT  and  t  are  written  in  terras  of  two  new  parameters,  cr 
3?  z  z  s 

and  ft ,  for  the  yield  condition,  a  and  ft  are  both  functions  of  r,  z,  t. 

s 


The  Flow  Conditions  for  the  Soil 


The  Equations  11-29,  11-30,  and  11-31  given  above  apply  to  the  so- 
called  "static"  rigid-plastic  yielding  case.  For  dynamic  loading  con¬ 
ditions  to  and  beyond  failure,  as  is  the  case  for  missile  penetration. 
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not  only  will  yielding  occur  but  also  viscous  flow  will  be  present. 
Recent  tests  have  indicated  that  the  angle  of  internal  friction  may  not 
be  influenced  much  by  the  rate  of  deformation,  but  the  cohesion  term 
appears  to  be  highly  influenced  by  the  deformation  rate. Por 
this  condition  the  Mohr-Coulomb  envelope  wouldbe  changed  and  expanded 
during  a  dynamic  flow  process  in  direct  proportion  to  the  change  in 
the  apparent  cohesion  term.  If  this  increase  in  apparent  cohesion  is 
assumed  to  be  proportional  to  the  rate  of  shear  deformation,  as  a  first 
approximation,  the  following  flow  condition  would  be  obtained. 


sf 


/av  av  \ 

=  C  +  K  ( -  1  ^  O'  'ton  0 


dz 


dr 


(11-32) 


where 


T  =  shear  resistance  on  the  flow  surface 

St  - 


C  =  cohesion  as  defined  before 


Kent  A.  Healy,  'U'ndrained  Strength  of  Saturated  Clayey  Silt, 

M.I.  T.  Report,  March  1963,  No.  15,  63-19,  Contract  No.  DA'^-079- 
Eng.  -224  with  U.  S.  Army  Engineers  Waterways  Experiment  Station, 
Dept,  of  Army. 

17 

A.  Casagrande  and  S.  D.  Wilson,  "Effect  of  Rate  of  Loading  on 
the  Strength  of  Clays  and  Shales  at  Constant  Water  Content,"  Geo- 
technique,  II  (1951),  pp.  251-263. 

1 8 

Robert  V.  Whitman  and  Kent  A.  Healy,  Shearing  Resistance  of 
Sands  During  Rapid  Loading.  M.  I.  T.  Report,  May  1962,  No.  9, 
62-113,  Contract  No.  DA-22 -79-Eng.  -224  with  U.  S.  Army  Engineers 
Waterway  Experiment  Station,  Dept,  of  Army. 
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a  =  effective  normal  stress  as  defined  before 
s 

</)  =  ultimate  angle  of  internal  friction  as  defined  before 

K  =  apparent  soil  viscosity  coefficient 

=  velocity  in  r  direction 

V  =  velocity  in  z  direction 
z 

a  V  a  V 

r  z 

- + - =  rate  of  shear  deformation 

az  ar 


If  C’  is  called  the  apparent  cohesion  then 

/av  av  X 

C  =  C  +  k( - -  + - -)  (11-33) 

\  az  ar  / 

The  coefficient  of  soil  viscosity  K  is  in  Equation  11-31  multiplied 

19 

by  the  rate  of  deformation  rather  than  the  velocity  gradient.  This  is 

justified  since  the  flow  condition  equation  is  actually  a  Bingham  body- 
20 

type  flow  formula  where  the  initial  yield  shear  resistance  is  a  function 
of  the  normal  stress  only.  A  plot  of  this  "law"  is  shown  in  Figure  10 


1 9 

Markus  Reiner^  Deformation,  Strain  and  Flow,  (London: 

H.  K.  Lewis  &  Co.  ,  Ltd.  ,  1960),  p.  20^ 

2  0 

A.  M.  Freudenthal,  The  Inelastic  Behavior  of  Engineering 
Materials  and  Structures,  (New  York:  John  Wiley  &  Sons,  1950), 
p.  402. 


Yield  limit 


Shear  resistance 


1 


^sy  =  ^  ^s  ^ 


+  cTg  tan  0  +  K 


Figure  10.  Shear  Resistance  vs  Rate  of  Deformation  for  a  Viscoplastic 
Material  with  a  Stress  Dependent  Yield  Limit 
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Substitution  of  the  flow  Equation  11-32  into  the  Equations  11-24 
through  11-26  obtained  from  the  Mohr  diagram  yields  a  set  of  equations 
similar  to  the  yield  Equations  11-29  through  11-31 


A  Mohr-Coulomb  plot  for  both  the  yield  and  flow  conditions  is  shown 
in  Figure  11. 

It  will  be  noticed  from  the  Mohr  diagram  that 
1.  T  is  positive  in  the  equation  of  motion  when 


TT 

4 


0  O  3  TT 


2 


T  is  negative  in  the  same  sense  when 


Sir 

~T 


1 

2 


2. 
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3.  CT  >  CT  when 
z  r 


4  2^2  2 


and  when 


6  5  TT  6 


4.  a  >  a  when 
r  z 


IT  6  ^  6 

-  -  I  <  13  <  TT  -  I 


For  special  angles  of  /S  the  following  conditions  hold 


TT  <f) 


/3  =T- 


4  2 


T  =0 
rz 


a  >  O’ 
z  r 


/3  =  — 

P  4  2 


T  =0 
rz 


a  >  a 
r  z 


/3  =  2 


a  >  (7  =  a 

z  r  s 


T  =  T 

rz  s 


j3  =  TT 


a  >  a  =  CT 
z  r  s 


-T  =  T 

rz  s 


P  = 


CT  >  CT  =  CT 

r  z  s 


T  3:  T 

rz  s 


j3  =  TT  -  (/) 


CT  >  CT  =  CT 

r  z  s 


-T  =  T 

rz  s 
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CT  =  a  T  =  T 

r  z  rz  max 

C7  =  CT  -T  =  T 

r  z  rz  max 

By  algebraic  manipulation,  formulas  can  be  obtained  for  a  non- 
dimensional  form  of  Equations  11-34  through  11-36.  Plots  of  these  non- 
dimensional  forms  are  shown  in  Figure  12,  using  the  argument  j3. 

The  evaluation  of  the  soil  viscosity  coefficient  during  dynamic 
loadings  must  be  done  experimentally.  It  is  a  basic  material  property. 
Research  now  is  being  conducted  at  the  University  of  Virginia  and  Sandia 
Corporation  to  evaluate  the  changes  in  (f>  and  C  during  dynamic  loadings 
using  a  simple  shear  device. 

The  equations  derived  for  the  shear  resistance  during  yielding 

and  flow  may  be  limited  under  certain  conditions,  because  it  is  also 

possible,  if  all  of  the  principal  stresses  are  high  and  applied  with  a 

very  high  rise  time,  to  break  the  grains  and  to  increase  the  surface 

21 

area  of  the  soil  particles.  This  is  called  comminution.  The  commi¬ 
nution  effect  is  seen  in  the  case  where  a  piston  is  shot  into  a  rigid  pipe 

23 

H.  L.  Morrison  and  W.  A.  Allen,  Projectile  Energy  E?q)ended 
in  Comminution  of  the  Target,  NOTS  1181,  NAVORD  Report  4885, 

U.  S.  Naval  Ordnance  Test  Station,  China  Lake,  California,  July  1955. 


'22 
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22 

filled  with  sand.  It  may  be  possible  in  such  a  test  not  only  to  unlock 
the  structure  of  the  soil,  then  break  the  grains?^  but  also  if  the  energy 
is  sufficient,  the  soil  might  even  be  liquified  or  vaporized.  The  energy 
required  to  comminute  soil  is  high.  During  initial  impact,  commi¬ 
nution  is  known  to  occur  in  the  velocity  regime  under  study^^  when 
the  target  is  composed  of  granular  materials.  Thus,  the  flow 
equations  developed  could  not  be  expected  to  hold  until  the  projectile 
has  truly  entered  the  soil  and  shear  flow  established.  The  time 
required  to  establish  shear  flow  is  a  function  of  the  velocity  of  wave 
propagation  in  soil. 

In  summary,  it  is  doubtful  that  the  yield  and  flow  equations 
would  hold  when  the  soil  were  being  comminuted,  liquified  or  vaporized. 

^  22 

M.  E.  Backman  and  H.  L.  Morrison,  The  Fracture  of 
Aggregates  of  Brittle  Particles,  NOTS  TM  No,  1582,  U.  S.~Naval 
Ordnance  Test  Station,  China  Lake,  California,  November  1953. 

23 

A.  J.  Hendron,  R.  E.  Fulton  and  Bijan  Mohraz,  The  Energy 
Absorption  Capacity  of  Granular  Materials  in  One -Dimensional 
Compression,  Air  Force  Special  Weapons  Center,  Kirtland  Air  Force 
Base,  New  Mexico,  January  1963. 

24 

Taylor,  ^  p.  354^ 
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Changes  in  Density  During  Shear 


As  soil  is  sheared,  it  expands  or  contracts  in  volume  so  that 

the  density  becomes  nearly  constant  depending  on  the  normal  stress 

perpendicular  to  the  surface  of  movement  and  stays  approximately  at 

25 

this  density  as  shear  distortion  is  increased.  This  ultimate  density 
is  called  the  critical  density.  Typical  plots  of  density  versus  defor¬ 
mation  are  shown  in  Figure  7b.  Critical  density  is  not  much  influenced 

by  the  state  of  stress  as  long  as  the  stresses  are  sufficient  to  cause 

26 

movement  along  yield  surfaces.  Arthur  Casagrande  was  the  first  to 

investigate  this  phenomenon  using  the  direct  shear  test;  later  Donald 

25  27 

Taylor  expanded  the  work  also  using  the  triaxial  test.  K.  H.  Roscoe 

has  continued  these  investigations  using  a  simple  shear  device. 

For  saturated  soils,  there  can  be  no  change  in  total  density, 
since  for  high  rates  of  shear  deformation  the  soil  does  not  have  time 
to  drain.  The  thwarted  changes  in  volume  of  the  soil  result  in  a 


Taylor,  loc.  cit. 
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Arthur  Casagrande,  "Characteristics  of  Cohesionless  Soils 
Affecting  the  Stability  of  Slopes  and  Earth  Fills,"  Contributions  to  Soil 
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corresponding  rapid  buildup  of  pore  water  pressure  and  a  corre¬ 
sponding  decrease  in  the  effective  normal  stress.  It  would  be  expected 
that  loose  saturated  soils  would  develop  high  positive  (compressive) 
pore  pressures  during  dynamic  loading  so  liquifaction  would  result; 
correspondingly,  dense  saturated  soils  would  develop  high  negative 
(tensile)  pore  pressures  during  shear  deformation  at  high  rates. 

Dry  and  partially  saturated  soils  are  usually  found  in  nature 
above  critical  density  or  below  critical  void  ratio.  An  exception  is 
loess  and  it  may  be  that  the  deformation  required  to  reach  critical 
density  is  higher  for  loess  than  for  other  soils  since  it  is  so  com¬ 
pressible,  Nevertheless,  some  critical  density  should  exist  for 
loess  after  its  structure  has  collapsed  by  deformation.  The  increase 
in  volume  for  dry,  dense  granular  materials  is  something  like  10  per- 
cent,  but  the  decrease  in  volume  of  unsaturated  loess  and  other  loose 
soils  may  be  somewhat  higher. 

If  the  initial  changes  in  density  are  ignored,  the  density  of  the 
soil  can  be  assumed  to  remain  constant  during  the  penetration  event, 
but  this  density  must  be  taken  as  the  critical  density  rather  than  the 
in  situ  density.  Then  in  the  continuity  equation,  the  density  could  be 
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considered  constant,  the  material  incompressible,  and  the  changes  in 
density  due  to  shear  deformation  can  be  ignored.  If  the  density  in  fact 
remains  constant,  it  would  be  expected  that  the  surface  of  soil  half 
space  around  the  point  of  impact  would  be  raised  and  that  the  volume 
of  the  heave  would  be  approximately  equal  to  the  hole  in  the  half  space 
caused  by  the  penetrating  missile. 

Relation  Between  the  Angle  of  Flow  and  the  Velocity  Copiponents 

In  the  flow  Equations  11-34  11-35  and  II- 36  t  ,  a  ,  and  o  are 
^  rz  r  z 

functions  of  cr  and  j3  at  any  instant  of  time.  C,  (j),  and  K  are  assumed 
s 

constant  for  a  given  soil  after  yield  conditions  are  reached.  The 
angle  /3  designates  the  direction  of  yield.  The  direction  of  maximum  a 
at  any  space  point  in  the  field  that  has  yielded  must  be  related  to  the 
velocity  components  by 

V 

tan  a  =  (11-37) 

z 

where  a  is  measured  counter  clockwise  from  the  z  axis.  V  and  V 

r  z 

are  the  components  of  the  maximum  velocity  in  the  r  and  z  directions. 

It  has  been  shown  that  the  direction  of  resultant  velocity  makes  an 
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angle  (j)  with  the  direction  of  the  yield  lines  because  plastic  deformation 

2  9 

must  be  accompanied  by  an  increase  in  volume  of  the  soil  if  0  >  0. 

30 

This  work  has  been  extended  so  as  to  calculate  a  velocity  field  for 

two  dimensional  plane  strain  problems  using  the  concept  of  a  plastic 
31 

potential.  But  these  velocity  fields  are  for  situations  where  the  soil 
reaches  a  critical  density  when  sheared  through  large  deformations. 
The  velocity  solutions  based  on  the  plastic  potential  would  probably 
hold  for  initial  motion  of  soil  particles;  but  with  continued  motion, 
it  will  be  assumed  that  the  direction  of  flown'  is  equal  to  j3  ,  the 
direction  of  yield.  Then  the  maximum  velocity  vector  is  in  the 
direction  of  active  slip,  and  the  maximum  shear  strain  rate  occurs 
in  the  direction  of  yield  or  active  slip. 


Thus 


V  =  V  sinS  +  V  cos  8 
s  r  z 


(11-38) 


tanjS  =  — 
z 


(11-3  9) 


4  a 
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where  V  =  the  magnitude  of  the  velocity  at  a  point  in  target 

s 

space. 


The  Radial  Normal  Stress  Assumption 


The  azimuthal  normal  stress  a  ,  has  some  influence  on  the  accel- 

eration  of  the  soil  as  shown  by  the  equation  of  motion  (II- 3).  The  o-^ 

normal  stress  however  does  not  influence  the  stresses  on  the  flow  or 

yield  surface  as  shown  by  the  calculations  on  page  33.  By  considering 

the  Boussinesq  equation  for  a  point  load  on  an  elastic  half  space,  it  is 

33 

seen  that  is  always  zero  or  tensile.  When  Poisson's  ratio  is  one- 
half,  as  for  an  incompressible  medixjm,  the  stress  cr^  is  zero.  But 

cannot  be  considered  equal  to  zero  since  soil  is  not  an  elastic  mate¬ 
rial.  The  maximum  and  minimum  principal  stresses  must  lie  in  the 
r-z  plane  for  the  axially  symmetrical  state  of  stress  since  movement 
can  only  be  in  the  r-z  plane.  Thus,  the  radial  normal  stress  must 

be  the  intermediate  principal  stress.  According  to  the  Haar-von 

34 

Karman  hypothesis,  the  limiting  state  of  stress  at  yield  is  reached 
when  two  of  the  three  principal  stresses — that  is,  when  two  Mohr 
circles — are  tangent  to  the  limiting  Mohr  envelope.  Then,  must  be 
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equal  to  either  the  maximum  principal  stress  or  the  minimum 
principal  stress.  If  a  radial  soil  element  is  being  moved  and  strained 
so  that  the  strain  in  the  direction  is  extension,  as  in  the  case  of 
movement  away  from  the  center  line  or  z  axis,  the  radial  normal 
stress  then  would  have  to  equal  the  minimum  principal  stress.  Then, 
from  Figure  9  it  can  be  seen  that  at  yield  or  flow 

ff  +  CT 

“■4,^  =  -  ’’s  (11-40) 

or  +  (C  +  tan*)  (11-41) 

Summary  of  Constitutive  Equations 

The  following  six  equations  can  now  be  added  to  the  equations 
from  continuum  mechanics.  The  boundary  conditions  and  initial  con¬ 
ditions  must  still  be  defined,  and  also  the  equation  of  the  surface  A 
of  the  projectile  and  the  mass  of  soil  carried  by  the  projectile  M 

must  be  determined,  since  two  new  parameters  a  and  ,8  have  been 

s 

introduced.  Henceforth,  a  ,  a  ,  a,,  a  ,  t  ,  and  t  will  be  con- 

r  z  4*  s  s  rz 

sidered  the  effective  stresses  at  flow  and/or  yield,  and  the  subscripts 
f  and  y  will  be  dropped  as  well  as  the  bar  to  denote  effective  stress. 


a 

z 


a  + 
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C  +  <7  tan  (^  +  K  I  ^ 
s  '  9z 


/9  V 
(■  ■' 
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2  cos  /3  sin(6  +  ^ 
cos  ^ 
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(11-34) 


(11-35) 


(11-36) 


(11-39) 


p  =  constant 


(11-42) 


CT  + 
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C  +  CT  tan  0  +  K  .  ^ 
s  \  8z 


/  9  V 

( 


+ 


9  V 
_ z 

9r 


)1 


sin  0  -  1 ' 
.  cos  0 


(11-41) 


If  CT  ,  CT  ,  CT,  and  T  are  less  than  the  right  hand  side  of  Equations  11-34, 
z  r  i//  rz 

n-35,  11-36,  and  11-41,  respectively,  then  the  material  is  considered 


to  be  rigid  and  no  movement  is  occurring. 


CHAPTER  III 


EXPERIMENTAL  INVESTIGATION 

I.  INTRODUCTION 

The  basic  purpose  of  the  experimental  program  was  (1)  to 
establish  the  boundary  and  initial  conditions  which  the  governing  dif¬ 
ferential  equations  must  satisfy  during  the  penetration  event;  (2)  to 
study  the  mechanism  of  how  the  soil  is  displaced  during  penetration 
so  that  realistic  assumptions  concerning  the  velocity  and  stress  fields 
at  any  instant  of  time  could  be  made;  and  (3)  to  obtain  evidence  of  the 
validity  of  the  theoretical  work. 

The  secondary  purposes  were  (1)  to  study  the  effect  of  projec¬ 
tile  dimensions,  weights,  and  impact  velocities  upon  the  mechanism  of 
penetration;  (2)  to  consider  the  effect  of  target  constraints  and  initial 
state  of  stress  in  targets;  (3)  to  obtain  some  acceleration  time  signa¬ 
tures,  using  controlled  laboratory  experiments,  that  would  give  some 
indication  of  the  gross  expenditure  of  the  projectile's  energy  at  any 
instant  of  time;  and  (4)  to  conduct  some  full  scale  tests  that  would  com¬ 
plement  the  laboratory  experiments  and  permit  correlation  with  labora¬ 
tory  tests. 
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II.  TYPE  OF  EXPERIMENTS 

There  were  two  kinds  of  laboratory  experiments — two  dimen¬ 
sional  and  three  dimensional.  A  two-dimensional  test  was  one  in  which 
the  target  material  was  confined  between  two  clear  Plexiglass  plates  so 
that  the  whole  event  of  penetration  could  be  recorded  using  high  speed 
photography.  The  thickness  of  the  target  medium  was  one  sixteenth  of 
an  inch  greater  than  the  thickness  of  the  projectile.  A  three-dimensional 
laboratory  test  was  one  in  which  the  right  circular  cylinder  was  shot 
into  a  tank  filled  with  the  target  medium. 

Gelatin  was  used  to  model  a  cohesive  soil  in  some  of  the  two- 
dimensional  tests.  Dry  Ottawa  sand  was  used  in  both  the  two-  and  three- 
dimensional  tests  to  model  a  noncohesive  soil.  Stacked  steel  needle 
bearing  rollers  were  also  used  in  some  two-dimensional  tests  to  model 
a  noncohesive  soil.  Water  was  used  in  some  three-dimensional  tests  as 
a  lower  bound  on  soils,  since  soils  are  composed  of  mixtures  of  soil 
particles,  water,  and  air. 

The  full  scale  tests  reported  were  conducted  at  (1)  Sandia  Cor¬ 
poration,  Area  III,  Albuquerque,  New  Mexico;  (2)  Sandia  Corporation, 
Tonopah  Test  Range,  Nevada;  and  (3)  at  Eglin  Air  Force  Base,  Florida. 
The  experiments  conducted  were  as  follows: 

1.  Two-dimensional  laboratory  tests: 

a.  Constant  velocity  tests  into  steel  needle  bearing  rollers 
using  different  nose  shapes,  projectile  widths  and 
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velocities.  Actually,  180  tests  were  ccnducted,  but 
only  a  summary  of  typical  tests  is  reported.  All  were 
tests  in  which  the  penetrator  was  moved  down  verti¬ 
cally. 

b.  Dynamic  penetration  tests  in  which  parallelepiped 
projectiles  were  shot  into  the  various  target  media. 

The  type  of  targets  and  number  of  tests  are  as  follows ; 

1.  Steel  needle  bearing — two  tests 

2.  Moist  bulked  silt — four  tests 

3.  Dry  sand — nineteen  tests 

4.  Gelatin — thirty-three  tests 

Results  from  only  one  of  each  of  the  above  tests  are  pre¬ 
sented  in  this  study.  Three  of  the  dry  sand  tests  and  four 
of  the  gelatin  tests  were  made  in  which  the  trajectory  of 
the  projectile  was  horizontal.  A  250-pound  pull  mechani¬ 
cal  actuator  was  used  in  the  gelatin  tests  as  a  velocity 
generator.  A  light  gas  gun  actuator  was  built  to  propel 
the  projectile  for  all  the  other  tests. 

2.  Three-dimensional  laboratory  tests: 

a.  Vertical  tests  into  water  using  a  right  circular  cylindri¬ 
cal  projectile:  High  speed  photography  was  used  to 
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record  the  motion  of  the  projectile  and  water.  Eleven 
tests  were  performed,  but  results  from  only  one  test  are 
presented  in  this  study. 

b.  Vertical  tests  into  sand  using  right  circular  cylindrical 

projectiles:  The  diameter,  weight  and  the  initial  velocity 
were  varied  as  well  as  the  density  of  the  sand.  One  hun 
dred  and  eleven  tests  were  performed,  but  results  from 
only  eleven  tests  are  reported.  Microwave  equipment 
was  used  to  record  the  motion  of  the  pi'ojectile. 

Several  other  tests  were  performed  using  various  nose 
shapes  at  high  velocities  generated  by  explosive-actuated 
devices.  No  motion  time  history  could  be  obtained,  but 
the  results  will  be  briefly  discussed. 

c.  Horizontal  tests  into  cylinders  of  pseudo  soils:  These 

tests  were  conducted  by  other  researchers  in  Division 
9226  at  Sandia  Corporation.  The  cylinders  of  pseudo 
soil  are  about  24  inches  in  diameter,  24  inches  long, 
and  composed  of  various  particles  mixed  with  a  plastic 
binder.  The  projectile  was  propelled  by  a  rocket  sled; 
projectile  nose, shape,  length,  and  velocity  were  varied. 
The  unconfined  compressive  strength  and  density  of  the 
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targets  were  varied.  One  hundred  and  twenty-one  tests 
have  been  conducted,  but  only  two  tests  will  be  dis  - 
cussed. 

3.  Pull  Scale  Tests : 

Essentially  vertical  tests  into  natural  deposits  of  soil: 
The  projectile  parametrs  (velocity,  weight,  diameter, 
length  and  nose  shape)  have  been  varied  and  many 
different  soils  have  been  penetrated.  Three  blunt- 
nosed  tests  into  a  dry  lake  playa  are  part  of  this 
study  and  will  be  reported.  Three  sharp-nosed  tests 
in  the  dry  lake  playa  and  one  into  clayey  sand  are 
reported,  but  these  tests  have  been  conducted  by 
others  at  Sandia  Corporation.  Time  vs.  deceleration 
data  are  reported  for  these  tests  when  the  accelerom¬ 
eters  and  telemetering  systems  survived. 
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III.  DESCRIPTION  OF  EXPERIMENTS 
Two-Dimensional  Laboratory  Tests 

a.  Constant  Velocity  Penetration  Tests  --  A  series  of 
continuous  vertical  penetrations  was  performed  using  a  blunt 
penetrator  at  various  constant  velocities  into  a  uniform  density 
granular  target  medium  having  no  cohesion. 

The  equipment  used  for  this  investigation  consisted  of  a  me¬ 
chanical  linear  actuator,  a  large  target  tank  with  a  transparent  face, 
an  array  of  steel  needle  bearing  rollers,  a  load  measuring  and  re¬ 
cording  system,  and  a  sequence  camera  for  recording  optical  data. 

See  Figure  13.  The  linear  actuator  used  was  driven  by  a  28-volt 
constant  speed  electric  motor.  A  pair  of  lead  screws  driven  by  the 
motor  through  a  gear  reduction  system  actuated  a  yoke.  The  pene¬ 
trator  was  mounted  on  the  yoke  and  was  guided  so  that  a  reproducible 
vertical  stroke  was  obtained. 

The  target  tank  was  46  inches  wide  by  36  inches  high.  It  was 
7 /8-inch  thick  so  that  the  3/4-inch  thick  penetrator  could  enter  without 
encountering  any  friction  on  its  front  or  back.  By  making  the  thick¬ 
nesses  (z  dimension)  of  both  the  target  material  and  the  penetrator  the 
same,  a  true  two-dimensional  (x  and  y)  test  was  achieved. 


Figure  13.  Setup  of  Two-Dimensional  Constant  Velocity  Penetration 
Test  Into  a  Simulated  Granular  Cohesionless  Medium 
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The  front  face  of  the  target  tank  was  1-inch  thick  Plexiglas  to 
permit  observation  of  the  movement  of  the  target  material  during 
penetration.  By  removing  this  front  face,  the  target  material  could 
be  restacked  and  returned  to  its  original  density,  and  a  grid  of  lines 
could  be  repainted  on  it  for  the  next  test.  Methods  were  developed  for 
these  operations  so  that  they  could  be  performed  quickly  and  with 
reproducible  results.  The  inside  surface  of  the  front  face  of  the  target 
tank  has  a  grid  of  black  lines  painted  over  half  of  it. 

The  target  material  consisted  of  many  thousands  of  precision- 
ground  steel  needle  bearing  rollers  of  three  diameters,  1/ 16 -inch, 

3/32 -inch  and  l/8-inch;  all  were  3/4-inch  long  and  had  rounded  ends. 
The  three  different  diameters  were  randomly  distributed  throughout 
the  array.  They  were  stacked  in  the  target  tank  so  that  the  length 
dimension  coincided  with  the  thickness  of  the  tank. 

To  study  movements  in  the  target  material  caused  by  penetration 
forces,  a  2 -inch  square  grid  of  narrow  white  lines  was  painted  on  the 
front  end  surfaces  of  the  stacked  rollers  before  each  test.  After  each 
test,  the  paint  was  carefully  removed,  the  array  of  rollers  was 
carefully  restacked  as  necessary  to  achieve  the  original  packing 
(density),  and  a  new  grid  was  painted  for  the  next  test. 
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The  load  measuring  and  recording  system  consisted  of  a  linear 
strain  gauge  connected  to  a  Brown  recorder.  The  strain  gauge 
assembly  was  calibrated  before  and  after  each  test  series.  The  Brown 
recorder  also  was  calibrated  immediately  before  and  after  each  test. 
The  recorder  operated  at  a  constant  paper  speed,  and  an  electric  clock 
reading  in  hundredths  of  a  second  was  included  in  the  pictures  taken  of 
the  target  tank  during  each  test. 

The  optical  observations  of  the  operation  of  each  test  were  made 
with  a  Nikon  Model  F  single  lens  reflex  35  millimeter  camera  equipped 
with  a  battery-powered  motor.  All  sequence  photographs  of  the  tests 
were  taken  at  the  fastest  frame  rate  of  the  camera,  between  three  and 
four  frames  per  second.  Kodak  Plus  X  black  and  white  film  was  used, 
and  full-frame  8-1/2  x  11 -inch  enlargements  were  printed  with  an 
automatic  photo-printer  at  a  constant  setting  to  provide  reproducible 
measurements  for  an  entire  individual  test  sequence. 

The  test  variables  investigated  included:  (1)  four  different  pene- 
trator  nose  shapes;  blunt,  triangular,  ogive,  and  semi-circular; 

(2)  three  different  penetrator  widths;  one  inch,  two  inches,  and  three 
inches;  and  (3)  five  different  penetrator  velocities:  1.  8,  2.  0,  2.  4,  3.  0, 
and  4.4  inches  per  second.  At  least  two  tests  were  conducted  using 
each  of  the  above  nose  shapes,  widths  and  velocities  for  a  total  of  one 
hundred  and  twenty  tests. 
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Figures  14,  15,  16,  and  17  show  a  series  of  sequence  photos 
taken  of  a  test  using  a  blunt  penetrator  three  inches  wide  at  a  constant 
penetration  velocity  of  2.  0  inches  per  second.  It  can  be  seen  that 
particles  greater  than  approximately  one  projectile  width  away  from 
the  projectile  centerline  moved  toward  the  surface  and  away  from  the 
projectile.  Particles  closer  to  the  projectile  than  this  were  first  dis¬ 
placed  toward  the  surface  due  to  the  flow  of  particles  from  beneath  the 
penetrator,  then  they  were  dragged  down  with  the  projectile,  but  they 
never  attained  a  depth  greater  than  their  original  depth.  The  distortion 
of  the  horizontal  grid  lines  indicated  that  a  boundary  layer  forms  along 
the  penetrator  in  this  velocity  range.  Further,  a  distinct  shear  front 
was  seen  to  move  through  the  medium  with  and  ahead  of  the  projectile 
in  much  the  same  manner  as  a  shock  wave  proceeds  in  front  of  a  high- 
velocity  projectile.  The  shear  front  is  defined  as  the  line  which  bounds 
the  zone  in  which  no  shearing  of  the  medium  has  occurred.  The  shape 
of  this  front  is  dependent  upon  projectile  velocity;  it  is  more  blunt  at 
lower  velocities.  The  location  of  the  shear  front  in  the  medium  caused 
by  the  penetration  event  can  clearly  be  seen  since  it  is  characterized  by 
a  narrow  band  of  material  having  an  increased  void  ratio.  On  the  sur¬ 
face  of  the  simulated  cohesionless  soil,  the  distance  from  the  center 
line  of  the  penetrator  to  the  shear  front  monotonically  increased  as  the 
depth  of  penetration  increased. 
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Typically,  particles  ahead  of  the  front  were  stationary.  Particle 
motion  started  as  the  front  passed  any  point.  If  a  given  grid  square  of 
particles  in  the  two-dimensional  test  is  examined,  the  square  is  seen 
to  distort  as  the  shear  front  passed  through  it.  Once  the  front  had 
completely  passed  through,  the  distorted  square  moved  as  a  rigid  body 
toward  the  surface  and  away  from  the  penetrator.  Horizontal  motion 
of  the  square,  however,  is  a  finite  bounded  quantity  for  any  penetration 
depth.  In  contrast,  the  vertical  displacement  increased  monotonically 
with  the  depth  of  projectile  penetration  because  the  penetrator  is  con¬ 
tinually  occupying  more  of  the  space,  pushing  the  particles  upward 
toward  the  direction  of  least  resistance. 

For  the  2-  and  3 -inch  wide  penetrate rs,  slip  surfaces  developed 
deep  in  the  medium.  These  are  characterized  by  stationary  material 
beneath  the  plane  and  the  material  above  the  slip  surface  moving  en¬ 
tirely  as  a  rigid  body.  As  noted  in  Figure  17,  the  grid  squares  above 
the  slip  surface  are  still  square  after  their  motion,  yet  they  have  been 
moved  up  and  away  from  the  penetrator  as  well  as  having  been  tilted. 

In  the  conventional  bearing  capacity  theory,  it  is  assumed  that 
a  wedge  of  the  penetrated  medium  forms  in  front  of  the  blunt  penetra¬ 
tor  and  stays  with  it.  This  wedge  is  seen  in  this  series  of  penetration 
tests,  but  the  particles  in  the  wedge  continually  change  with  increasing 
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depth  of  penetration.  A  wedge  of  the  same  particles  still  may  exist 
if  the  medium  has  cohesive  strength. 

A  preliminary  analysis  of  unit  resistance  at  a  depth  of  15.4  inches, 
using  three  velocities,  three  widths  and  four  nose  shapes  is  shown  in 
Figure  18.  The  curves  show  the  variation  in  average  unit  resistance 
(112  tests)  versus  the  length  to  diameter  of  the  four  nose  shapes. 

The  nose  shapes  tested  were,  blunt,  semicircular,  triangular  and  six 
caliber  ogive.  The  relationship  between  resistance  and  depth  of  pen¬ 
etration  was  nearly  linear. 

The  unit  resistance  to  penetration  decreases  with  increasing  pen- 
etrator  width  for  each  nose  shape  tested.  Intuitively,  this  situation  does 
not  seem  reasonable;  however,  Hvorslev  and  Kerisel  report  similar 
observations  under  considerably  different  test  conditions  but  at  approxi¬ 
mately  the  same  penetration  velocities.  Earlier  tests  in  the  same 
37 

medium  by  Colp  showed  this  same  result  for  different  nose  shapes. 
Unit  penetration  resistance  decreased  as  the  ^  ^  j  ratio  increased. 

— 

M.  J.  Hvorslev,  Discussion  at  Session  II  -  Shallow  Foundations, 
Symposium  on  Bearing  Capacity  and  Settlement  of  Foundations,  Duke 
University,  April  1965. 

Of' 

J.  Kerisel,  "Deep  Foundations  Basic  Experimental  Facts," 

Deep  Foundations  Conference,  Mexico  D.  F. ,  Mexico,  December  1964. 

07 

J.  L.  Colp,  An  Experimental  Investigation  ^  t^  Continuous 
Penetration  of  a  BlunTSody  into  a  Simulated  Cohesionless  Soil,  Sandia 
Corporation,  Albuquerque,  New  Mexico,  SC-RR-65-260,  December  1965. 


AT  DEPTH  15.4 

AVERAGE  UNIT  PENETRATOR  LOAD  (LB/ IN) 


76 


b.  Dynamic  Penetration  Tests  — 

( 1 )  Steel  needle  bearing  roller  penetration  tests:  Two 
dynamic  penetration  tests  into  a  two-dimensional  target  array  of  steel 
needle  bearing  rollers  were  performed.  These  tests  were  performed 
to  correlate  the  observations  made  in  the  low,  constant  velocity  pene¬ 
trations  described  above  to  high  velocity  penetrations.  In  this  series, 
the  velocity  generator  used  was  a  light  gas  gun  and  the  target  tank 
and  medium  were  the  same  that  were  used  in  the  low  constant  velocity 
series.  Motion  pictures  of  the  penetration  were  taken  at  a  rate  of 
1000  frames  per  second.  Figure  19  shows  the  test  set-up  used. 

The  first  shot  had  a  projectile  impact  velocity  of  225  fps. 

A  severe  shear  zone  developed,  in  addition  to  a  slip  surface,  well 
beneath  the  projectile.  The  medium  above  the  surface  moved  as  a  rigid 
body  and  maintained  the  respective  orientation  of  the  painted  grid 
lines.  However,  as  the  crater  formed,  the  slip  plane  became  a 
barrier  beyond  which  no  particles  of  the  medium  were  thrown  out  by 
the  impact.  Thus,  it  could  be  assumed  that  the  locking  bonds  of  the 
medium  which  moved  as  a  rigid  body  along  the  slip  plane  had  been 
weakened,  while  those  particles  beneath  the  plane  remained  strong 
and  provided  a  lower  limit  to  crater  size. 


Illi 


3t  Into  Steel  Roller  Target 
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The  photographic  coverage  of  the  shot  conducted  at  350  fps 
was  not  as  good  as  the  225  fps  shot  coverage,  and  the  existence  of  a 
defined  slip  surface  could  not  be  determined. 

The  projectile  penetrated  about  5  inches  into  the  roller 
array  for  both  impact  velocities.  No  scratching  of  the  Plexiglas  con¬ 
tainer  walls  could  be  seen,  so  it  can  be  assumed  that  the  entire  decel¬ 
eration  effort  was  provided  by  the  medium  itself.  None  of  the  rollers 
in  the  area  of  the  container  walls  were  disturbed,  so  it  is  not  believed 
that  the  container  afforded  false  boundary  conditions.  The  strength  of 
the  deceleration  forces  is  demonstrated  by  the  severe  plastic  deforma¬ 
tion  of  the  steel  projectile  shown  in  Figure  20.  At  the  time  of  impact, 
the  projectile  started  to  tumble  because  of  its  short  length  and  because 
it  did  not  impact  exactly  perpendicularly.  The  most  intense  deceleration 
occurred  within  the  first  half  inch  of  penetration.  It  is  quite  evident 
that  the  entry  phase  deserves  considerable  attention  in  the  penetration 
problem  analysis. 

In  the  low  constant  velocity  penetrations,  the  penetration 
force  during  the  first  half  inch  of  penetration  was  very  low.  For  this 
reason,  it  is  not  believed  that  resistance  force  correlation  between  low 
velocity  and  high  velocity  shots  can  be  made. 
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Mild  SLccJ 
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Steel  Roller  Indentations 
on  Projectile  Face  After 
350  ft/sec  Impact 


Sample  Steel  Rollers 
from  Target  Medium 


Figure  20.  Light  Gas  Gun  Projectile  After  Impact  Into 
Two-Dimensional  Steel  Roller  Target 
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Figures  21  through  28  show  individual  frames  from  the 
high  speed  movies  taken  of  this  test  series.  The  displacements  and 
motions  of  the  target  material  particles  during  and  after  the  projectile 
penetration  can  be  seen. 

The  velocity  of  the  particles  within  the  target  medium  dur¬ 
ing  the  penetration  event  is  quite  small;  in  most  cases  examined,  it 
was  about  one -fifth  the  magnitude  of  the  penetrator  velocity.  In  two 
high  velocity  penetrations  (350  and  225  fps)  into  the  roller  bearings, 
relatively  few  particles  in  the  general  area  of  the  penetration  event 
were  affected  until  after  the  projectile  stopped.  At  that  time,  rollers 
were  thrown  out  of  the  medium  at  a  velocity  about  one -fifth  that  of  the 
impacting  projectile  velocity.  Only  a  few  individual  rollers  could  be 
followed  in  the  motion  pictures  of  the  event,  but  their  exit  velocities 
appeared  to  be  quite  similar  in  magnitude. 

Further,  in  this  test  series,  the  particles  were  not  raised 
a  large  distance.  By  far  the  greatest  amount  of  energy  available  for 
the  penetration  process  is  required  to  break  the  bonds  of  the  locking 
granular  material  and  to  continue  the  shearing  process.  For  this 
reason,  it  is  proposed  that  in  any  preliminary  theory  of  medium  failure 
during  a  penetration  event,  only  the  shearing  action  need  be  considered. 
The  greater  sophistication  resulting  from  inclusion  of  terms  to  account 
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re  21.  High  Speed  Movie  Frame  of  Dynamic  Penetration  Into 
Steel  Roller  Target  -  Time  0.  0006  Second  After  Impact 


Figure  25.  High  Speed  Movie  Frame  of  Dynamic  Penetration  Into 
Steel  Roller  Target  -  Time  0.  0196  Second  After  Impact 
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Figure  28.  High  Speed  Movie  Frame  of  Dynamic  Penetration  Into 
Steel  Roller  Target  -  Time  0.  3066  Second  After  Impact 
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for  particle  exit  velocities,  etc. ,  should  be  added  only  after  a  much 
better  grasp  of  the  basic  shearing  action  is  attained. 

(2)  Moist,  bulked  silt  penetration  tests :  Four  low  velocity 
penetration  tests  were  made  into  a  moist,  bulked  silt  target.  A  direct 
acting  impact  mechanical  actuator  was  used  on  this  series.  The  target 
tank  used  was  similar  to  that  described  for  the  constant  velocity  steel 
roller  penetration  tests.  The  moist  silt  target  material  was  placed  in 
the  target  tank  in  a  bulked  state.  The  density  of  the  target  material  was 
uniform  throughout  the  target  tank.  The  target  was  constructed  by  plac¬ 
ing  measured  weights  of  material  in  measured  volumes  between  the 
Plexiglas  plate. 

High  speed  movies  were  taken  of  the  penetration  event  into 
the  bulked  silt.  Figures  29,  30,  and  31  are  prints  of  individual  frames 
from  the  high  speed  movies  taken  during  this  series.  The  decided  dif¬ 
ference  in  the  target  material  deformation  pattern  observed  in  similar 
tests  into  dense,  granular  materials  can  be  clearly  seen  in  these  photo¬ 
graphs. 

(3)  Dry  Ottawa  sand  penetration  tests :  A  series  of  pene¬ 
tration  tests  were  performed  into  a  target  of  clean,  dry  Ottawa  sand 
that  passed  a  20-mesh  screen  and  was  retained  on  a  30-mesh  screen. 
Nineteen  tests  were  included  in  this  series. 


Figure  29.  Initial  High  Speed  Movie  Frame  of  Dynamic  Penetration 
Into  Bulked  Silt  Target 


Figure  31.  Final  High  Speed  Movie  Frame  of  Dynamic  Penetration 
Into  Bulked  Silt  Target 
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A  light  gas  actuator  was  used  to  propel  the  projectile  for 
this  series.  The  target  tank  was  the  same  as  described  earlier.  Fig¬ 
ure  32  shows  the  test  setup  used. 

The  test  series  included  projectile  impact  velocities  from 
106  feet  per  second  to  478  feet  per  second.  Two  sizes  of  projectiles 
were  used;3/4*  x  3/4"  in  cross-section  and  1-1/2"  x  1-1/2"  in  cross- 
section.  All  projectiles  weighed  0.  51  pound.  The  density  of  the  Ottawa 
sand  target  material  was  uniform  for  each  individual  test  and  was  varied 
between  105  and  112  pounds  per  cubic  foot. 

High  speed  movies  were  taken  of  each  test  run.  Figures  33 
and  34  show  individual  frames  from  the  movies  taken.  As  can  be  seen, 
the  deformation  patterns  observed  are  quite  similar  to  those  shown  in 
the  constant  velocity  penetration  into  steel  needle  bearing  rollers  shown 
earlier. 

(4)  Gelatin  penetration  tests:  A  series  of  penetration  tests 
into  a  target  of  gelatin  were  performed.  Gelatin  was  selected  as  a 
target  material  to  simulate  a  soft,  viscous,  incompressible  cohesive 
material  such  as  saturated  clay.  Gelatin  is  a  double -refracting  ma¬ 
terial.  A  study  was  begun  in  which  photoelastic  means  were  used  to 
determine  the  internal  stresses  induced  by  penetration.  A  large 
polariscope  (Figure  35)  was  constructed  and  used  in  this  study.  Thirty- 
three  penetration  tests  into  this  target  material  were  made. 


Figure  32.  Laboratory  Test  Setup  for  Dynamic  Penetration 
Into  Ottawa  Sand 
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Figure  35.  Laboratory  Test  Setup  for  Penetration  Test  Into  Gelatin 
Showing  Large  Polaris  cope  Used 
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A  linear  mechanical  actuator  was  used  to  propel  the  pro¬ 
jectiles.  The  same  target  tank  as  used  in  the  previously  described 
needle  bearing  roller  tests  was  used  for  this  series.  Projectile  im¬ 
pact  velocities  were  varied  from  110  feet  per  second  to  395  feet  per 
second.  The  same  two  projectile  sizes  as  used  in  the  Ottawa  sand 
tests  described  above  were  used  in  this  series.  The  gelatin  was  of 
uniform  density  for  each  individual  test.  Due  to  the  difficulty  of 
maintaining  uniform  laboratory  temperature,  it  was  not  possible  to 
hold  the  gelatin  consistency  exactly  constant  from  test  to  test. 

High  speed  movies  of  the  penetration  into  the  gelatin  were 
taken  through  the  large  polariscope  used  in  this  test  series.  Figures 
36  through  38  show  individual  frames  from  one  of  the  high  speed  movies 
taken.  These  photos  show  the  double -refracting  effects  in  the  gelatin 
caused  by  the  entrance  of  the  projectile.  The  resulting  pictures  were 
of  such  poor  quality  that  the  stress  field  could  not  be  analyzed.  The 
principal  problem  was  that  not  enough  light  could  be  transmitted  through 
the  polariscope  in  the  short  exposure  time  required  for  1000  frames  per 
second  movies.  Three  different  lighting  systems  were  tried  and 
camera  speeds  up  to  25,  000  frames  per  second  were  employed.  It  is 
hoped  that  this  work  can  be  extended  in  the  future  using  a  series  of  ex¬ 
ploding  wires  as  light  sources.  Nevertheless,  the  mechanism  of  failure 
could  be  observed  with  these  pictures.  They  show  that  there  was  never 
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any  side  friction  on  the  projectile  during  the  event  except  when  the 
projectile  tumbled.  The  stresses  were  low  except  near  the  front 
of  the  projectile  and  along  the  sides  of  the  projectile.  The  initial  shock 
wave  in  the  target  outran  the  projectile.  Once  the  gelatin  yielded 
at  the  surface,  the  energy  going  into  the  shock  wave  was  terminated. 

The  direction  of  yielding  and  flow  at  any  instant  of  time  was  similar  to 
the  sand  and  needle-bearing  tests. 

Three-Dimensional  Laboratory  Tests 

a.  Water  Tests  --  A  series  of  11  experiments  was  conducted 
using  water  as  the  target  medium.  It  was  reasoned  that  the  behavior 
of  some  soils  in  their  saturated  states  might  approach  that  of  water. 
Most  natural  deposits  of  soil  are  mixtures  of  soil  particles,  air,  and 
water.  Saturated  clays  of  low  density  could  possibly  be  characterized 
as  a  viscous  fluid  during  shear  flow;  thus,  water  penetration  should 
give  a  lower  bound  on  the  resistance  to  penetration  for  incompressible 
soils. 

The  model  projectile  used  in  these  tests  was  1-1/2  inches  in 
diameter,  2-1/4  inches  long  and  weighed  l/2  pound.  It  was  shot  down¬ 
ward  into  a  glass  tank  of  water  and  high  speed  photography  was  used  to 
record  the  motion.  In  one  case,  the  projectile  was  allowed  to  fall  free. 
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The  velocity  generator  was  the  light  gas  gun  used  for  the  two- 
dimensional  tests.  The  water  tank  was  about  2  feet  square  and  about 
six  feet  deep.  The  entrance  velocities  used  for  these  11  tests  varied 
between  16  and  218  feet  per  second.  The  position  of  the  projectile  was 
plotted  versus  time  by  using  the  timing  marks  on  the  high  speed  photo¬ 
graph  frames.  By  numerically  differentiating  the  motion  vs.  time  plot, 
velocity  was  calculated.  By  differentiating  velocity  vs.  time,  accelera¬ 
tion  was  calculated.  A  typical  set  of  data  is  shown  for  test  number  2 
in  Figure  39.  A  picture  of  the  projectile  during  water  penetration  is 
shown  in  Figure  40.  In  all  cases,  there  was  a  cavity  around  the  pro¬ 
jectile  sides  and  rear,  and  this  cavity  stayed  with  the  projectile  until 
it  came  to  rest  on  the  bottom  of  the  tank.  Steady-state  hydrodynamic 
theory  assumes  the  resisting  force  to  movement  is  proportional  to  the 
square  of  the  velocity,  or  that  the  velocity  squared  divided  by  the 
deceleration  is  equal  to  a  constant  for  a  given  projectile.  In  all  water 
tests,  the  average  deceleration  was  not  proportional  to  the  entrance 
velocities.  Also  for  any  individual  test,  the  deceleration  was  not 
proportional  to  the  velocity  squared  at  any  instant  of  time.  In  the  data 
shown  for  test  number  2,  the  velocity  squared  divided  by  the  decel¬ 
eration  increased  36  percent  in  the  10  milliseconds  it  took  the  pro¬ 
jectile  to  travel  12-1/4  inches.  This  is  typical  of  all  of  the  tests. 


DIAMETER  OF  PROJECTILE - 1.5  INCHES 
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b.  Three-Dimensional  Sand  Tests  --  A  series  of  111  tests  was 
conducted.  The  tests  were  similar  to  the  water  tests  except  dry 
Ottawa  sand  at  various  densities  was  used  as  the  target.  Of  course, 
the  cameras  could  not  be  used  to  record  the  motion.  At  first  an 
attempt  was  made  to  measure  the  deceleration  of  the  projectile  directly. 
Accelerometers  were  mounted  in  the  projectile,  but  in  only  two  cases 
was  it  possible  to  record  a  part  of  the  deceleration  time  signature 
since  the  hard  wire  connection  to  the  recording  unit  invariably  broke. 

At  the  time  the  tests  were  being  conducted  there  existed  no  radio 
telemetering  package  small  enough  to  fit  into  the  1-1/2  inch  diameter 
projectiles,  so  a  microwave  system  was  developed  for  this  study  by 
Dr.  Don  Thorn  of  the  University  of  New  Mexico  to  measure  displace¬ 
ment  versus  time  during  the  event.  An  antenna  was  placed  in  the  bottom 
of  the  soil  container.  Microwaves  were  projected  upward  through  the 
soil  and  reflected  off  of  the  projectile  during  penetration.  There  was  a 
phase  difference  between  the  projected  and  reflected  signal.  The 
recording  system  allowed  the  superimposing  of  a  time  signal  on  the 
microwave  record.  The  emitted  and  reflected  signals  were  super¬ 
imposed  causing  either  cancellation  or  reinforcement.  The  time  be¬ 
tween  peaks  of  reinforcement  of  signal  was  measured;  then  using  the 
precalibrated  velocity  of  the  microwaves  through  the  sand,  the  distance 
moved  during  the  time  interval  was  calculated.  No  change  of  velocity 
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of  the  microwaves  could  be  detected  when  the  density  of  the  soil  was 
varied  between  99  and  110  pounds  per  cubic  foot.  When  the  projectile 
stopped,  the  signal  became  constant.  The  distance  penetrated  was 
measured.  The  initial  velocity  was  calculated  independently,  by 
measuring  the  time  increment  that  elapsed  for  the  projectile  to  beak 
two  carbon  leads  spaced  a  known  distance  apart  across  the  flight  path 
of  the  projectile  before  impact.  The  time  increment  was  measured 
with  a  Berkley  counter.  Working  back  from  the  maximum  depth  of 
penetration,  it  was  possible  to  calculate  the  entrance  velocity  and  total 
time  of  trajectory  in  the  soil.  In  only  11  of  the  111  tests  was  this 
possible  since  usually  the  projectile  tumbled  during  the  trajectory. 

At  first  the  distance  time  data  were  numerically  differentiated  to 
obtain  the  velocity  and  deceleration  plot,  but  this  resulted  in  rough 
erratic  plots  so  a  polynomial  was  fitted  to  the  depth-time  data  using  a 
computerized  method  to  best  fit  the  data.  These  functions  were  dif¬ 
ferentiated  twice  to  obtain  velocity  and  deceleration  plots.  These  data 
for  the  11  tests  are  shown  in  Figures  41  through  51.  The  4  tests  in 
which  the  unit  weight  of  the  soil  was  107  Ib/cu.  ft.  are  the  best  tests. 

In  none  of  these  cases  was  the  deceleration  proportional  to  the  velocity 
squared.  Because  the  projectile  tumbled,  there  is  no  uniform  trend  in 
the  relationships  between  deceleration  and  depth. 


DIAMETER  OF  PROJECTILE  —  0.75  INCHES 


108 


N0liVH3T303a 


DIAMETER  OF  PROJECTILE  —  0.75  INCHES 
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Figure  43.  Motion -Time  and  Motion -Penetration  Relations  for  Test  No.  (SO)G-IV-12A 
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UNIT  WEIGHT  OF  SAND—  lOT  POUNDS  PER  CUBIC  FOOT 
DIAMETER  OF  PROJECTILE  —  0.75  INCHES 
WEIGHT  OF  PROJECTILE — 0.5  POUNDS 
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UNIT  WEIGHT  OF  SAND— 107  POUNDS  PER  CUBIC  FOOT 
DIAMETER  OF  PROJECTILE—  1.5  INCHES 
WEIGHT  OF  PROJECTILE— 0.5  POUNDS 


114 


CD 

3 

£ 

<n  u 

1  i  I 

2  «q  2 


s 

£ 


X 

o 

u 

» 


^  UJ 


£ 

2 

a 


£ 

2 

0. 


K  O 
UJ 

hi  X 

Z  S 

<  UJ 

o  > 


n  u 

UJ 

(A 


O 


< 

00 

CO 


ii 


CQ 

<U 

H 

«w 

o 


m 

a 

o 

•rH 

4-» 

cd 

0) 

K 

c 

o 


cd 

!h 

■<-> 

0) 

a 

0) 

D. 


c 

o 


TO 

c 

oS 

(U 

a 

•iH 

H 

I 

c 

o 


t- 

<u 

u 

3 

bc 

•pH 


DIAMETER  OF  PROJECTILE — 1.5  INCHES 


115 


,03S/ld  5_OI 


Nouvudidoao 


DIAMETER  OF  PROJECTILE - 0.75  INCHES 


116 


o 

o 


u 

S 

3 

U 


o 

z 

3 

o 

a. 


u 

UJ 

o 

o: 

a. 


z 

o 

ui 


inaj  Nl  N0llV8i3N3d 


00 


o 

CD 

<u 

H 

Li 

O 

tw 

CD 

C 

o 


K 

C 

o 

•IH 

CO 

L 

-►» 

0) 

c 

0) 

Oh 

I 

c 

o 


TJ 

c 

cO 


H 

I 

c 

o 

•rH 

o 


02 

cu 

Li 

P 

M 

*rH 


119 


Soil  noses  were  found  on  the  projectile  after  the  tests,  and  the 
other  comminuted  soil  was  found  in  the  general  area  where  the  ejecta 
fell  back  in  the  tank.  The  soil  was  warm  to  the  touch  immediately 
after  each  test. 

Later,  three  different  1-1/2 -inch  diameter  projectiles  were  shot 
into  Ottawa  sand  using  an  explosive -actuated  velocity  generator  and 
longer  projectiles.  This  system  gave  a  stable  trajectory  and  velocities 
up  to  1100  feet  per  second.  Three  different  projectiles  were  used. 

The  shapes  of  the  noses  of  the  projectiles  were  blunt,  hemispherical 
and  2.  2  caliber  ogive.  In  these  tests  and  previous  tests,  a  nose  of 
pulverized  sand  was  found  on  the  blunt  nosed  projectile  and  the  hemi¬ 
spherical  nose  projectile,  and  a  partial  nose  was  found  on  the  sharp 
nosed  projectile.  Photographs  of  the  projectiles  are  shown  in  Figure 
52.  Photographs  of  the  soil  noses  are  shown  in  Figures  53  through  55. 
The  sides  of  the  projectile  were  still  smooth  after  the  event.  The  blunt 
nosed  hard  steel  projectile  was  very  rough  and  pitted  on  the  flat  face 
after  the  event.  The  sand  was  hot  to  touch  after  the  projectile  had  been 
fired  into  it. 

c.  Horizontal  Tests  Into  Pseudo  Soils  --  A  series  of  121  tests 
was  conducted  by  Division  9226,  Sandia  Corporation,  and  the  results 
are  discussed  briefly  here.  A  radio  telemetry  system  3/4  inch  in 
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a.  1-1/2  inch  blunt  nose  taken  off  cylindrical  afterbody  of 
projectile  shot  into  3-D  Ottawa  sand  target. 

b.  3/4  X  3/4  inch  parallelepiped  projectile  shot  into  2-D 
roller  bearing  target. 

c.  1-1/2  inch  cylindrical  projectile  with  hemispherical 
nose  shot  into  3-D  Ottawa  sand  target. 

d.  I-I/2  inch  cylindrical  projectile  with  six  caliber 
tangent  ogive  nose  shot  into  3-D  Ottawa  sand  target. 


Note: - line  of  separation.  See  Figures  53  through  55 

for  soil  noses  that  adhered  to  these  projectiles. 


Figure  52.  Photograph  of  Projectiles  After  Impact 
Showing  Nose  Distress 


DIAMETER  OF  PROJECTILE  —  1.5  INCHES 
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Figure  51.  Motion-Time  and  Motion -Penetration  Relations  for  Test  No.  (SO)G-IV-3D 


Figure  54.  Three  Different  Views  of  Soil  Nose  that  Formed 
on  Hemispherical  Nose  of  1-1/2  inch  Projectile. 

Ottawa  Sand  Target 
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diameter  and  about  4  inches  long  was  developed  which  included  a 
piezoresistive  accelerometer.  The  decelerations  were  transmitted 
during  penetration  and  recorded  on  magnetic  tape.  The  instrumenta¬ 
tion  seemed  to  be  excellent  and  should  in  the  future  provide  valuable 
data.  The  entrance  velocities  were  recorded  using  high  speed  cameras. 
The  procedure  was  to  propel  the  projectile  down  a  tube 
using  J.  A.  T.  O.  rocket  power.  The  propelling  equipment  was  stopped 
about  6  feet  short  of  impact  to  allow  the  projectile  to  fly  free  into  the 
target.  The  target  was  a  mixture  of  various  particles  and  plastic  en¬ 
closed  in  a  "sono"  tube.  The  diameter  of  the  tube  was  24  inches  and 
the  length  about  24  inches.  Two  sets  of  data  are  shown  in  Figures  56 
and  57.  It  has  been  noticed  that  the  deceleration  records  are  highly 
influenced  by  the  angle  of  attack.  If  the  projectile  was  not  exactly  per¬ 
pendicular  to  the  target,  the  entrance  spikes  were  lower  and  the  depth 
of  penetration  greater  for  the  blunt  nosed  projectiles.  The  sharp  nosed 
projectiles  were  hot  when  taken  from  the  target  but  the  blunt  nosed  pro¬ 
jectiles  were  cool.  Soil  noses  were  formed  on  the  blunt  projectiles.  In 
Figures  58  and  59  are  photographs  of  l-l/2-inch  blunt  nosed  projectiles 
flying  out  of  the  rear  of  the  targets  with  their  soil  noses  still  attached. 
The  density  of  the  target  material  was  56  pounds  per  cubic  foot  and  the 
unconfined  compressive  strength  was  3  pounds  per  square  inch. 
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Figure  58.  Blunt  Nosed  Projectile  Flying  Out  Rear  of  Horizontal  Target 
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The  data  obtained  so  far  generally  show  that  in  these  materials 
the  deceleration  is  almost  constant  as  the  projectile  passes  through 
the  constrained  target.  There  is  not  much  difference  in  these  decelera¬ 
tions  for  either  the  blunt  or  sharp  nosed  projectiles.  However,  the 
sharp  nosed  projectiles  do  not  show  the  high  entrance  spike  the  blunt 
nosed  projectiles  exhibit.  The  targets  after  penetration  have  radial 
cracks,  a  small  crater  at  the  entrance  point  and  a  much  larger  crater 
at  the  point  of  exit. 

Full  Scale  Tests 

a.  Blunt  Nosed  Tests  Into  Dry  Lake  Playa  --  A  series  of  three 
blunt  nosed  vertical  tests  has  been  performed  as  a  part  of  this  study. 

The  projectiles  were  dropped  from  a  helicopter  into  a  dry  lake  playa 
at  the  Sandia  Test  Range,  Tonopah,  Nevada.  In  Table  I,  the  pertinent 
projectile  parameters  and  depths  of  penetration  are  given.  In  Test 
No.  120-57  the  accelerometers  and  radio  transmission  package  inside 
the  projectile  were  shattered  upon  impact  so  no  deceleration  time 
signature  was  obtained.  In  the  two  other  blunt  nosed  tests,  a  record 
was  obtained  using  a  shock  mitigated  package,  but  the  initial  decel¬ 
eration  spike  was  not  obtained.  The  average  deceleration  after  the 
spike  was  approximately  2000  g's.  In  each  of  the  three  tests,  the 
heave  of  the  surface  was  only  about  l/4  inch  except  near  the  edge  of 
the  impact  crater,  but  the  volume  of  heave  plus  the  volume  of  the  crater 
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above  the  initial  surface  was  a  little  greater  than  the  volume  of  the 
hole  created  by  the  projectile. 

The  surface  elevations  after  penetration  are  given  in  Table  II 
for  Test  No.  120-57.  The  surface  deformations  for  Test  Nos.  120-89 
and  120-97  were  similar.  Table  III  shows  the  calculations  of  the 
volume  changes  for  Test  No.  120-57.  The  crater  was  larger  in  Test 
Nos.  120-89  and  120-97  than  in  Test  No.  120-57.  The  soil  strata 
descriptions  are  shown  for  Test  No.  120-57  in  Figure  60  and  for  the 
other  two  tests  in  Figure  61. 

The  ground  surface  around  the  point  of  impact  was  cracked 
radially  out  about  12  feet  in  Test  No.  120-57  and  somewhat  less  in  the 
other  two  tests.  In  Figures  62  and  63  are  shown  the  surfaces  after 
impact  before  and  after  the  cracks  were  marked  (Test  No.  120-57). 

A  soil  nose  formed  in  all  three  tests.  In  Figures  64a,  b,  and  c 
are  shown  the  soil  nose  taken  from  the  projectile  in  Test  No.  120-57. 

In  Figures  65a  and  65b  are  shown  the  projectile  before  removal  and  the 
soil  nose  being  excavated  after  the  projectile  was  removed.  In 
Figures  66a  and  66b  are  shown  photographs  of  the  projectiles  used 
in  Test  Nos.  120-87  and  120-97  after  removal  from  the  soil.  The 
paint  was  slightly  removed  on  the  side  of  the  projectile  that  was  on 
the  bottom  side  of  the  trajectory  but  this  removal  area  does  not  extend 
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TABLE  III 


Calculations  for  Volumes  of  Movement 
Test  No.  120-57 


Seg 

No. 

w 

w 

7 

w 

3 

^2 

^l"*'^2 

hj^+2h2 

h2^+2h2 

hj^+h2 

A^^Chj^+hj) 

y 

R  =  3cfy 

Vol  =  -jRA 

North 

1 

2 

3 

4 

5 

6 

0.137 

1.0 

1.4 

1.6 

4.8 

2.65 

0.07 

0.5 

0.7 

0.8 

2.4 

1.32 

0.05 

0.33 

0.47 

0.53 

1.6 

0.88 

0 

0.08 

0.03 

0.02 

0,01 

0.01 

0.08 

0.03 

0.02 

0.01 

0.01 

0 

0.08 

0.11 

0.05 

0.03 

0.02 

0.01 

0.16 

0.14 

0.07 

0.04 

0.03 

0.01 

2.0 

1.27 

1.40 

1.33 

1.50 

1.0 

0.005 

0.055 

0.035 

0.024 

0.048 

0.013 

0.091 

0.424 

0.658 

0.705 

2.400 

0.883 

1.06 

1.2 

2.2 

3.6 

5.2 

10.0 

1.154 

1.624 

2.858 

4.305 

7.600 

10.883 

0.01 

0.14 

0.16 

0.16 

0.573 

0.226 

Total  Volume  of  Heave  in  North  Quadrant 

(cu  ft) 

1.268 

Crater 

0.69 

0.34 

0.23 

0.42 

0 

0.42 

0.42 

1.0 

0.143 

0.23 

0.375 

0.605 

0,136 

South 

1 

0.47 

1.00 

1.09 

3.13 

2.98 

2.57 

11.24 

0.24 

0.50 

0.55 

1.56 

1.49 

1.28 

5.62 

0.16 

0,33 

0.36 

1.04 

0.99 

0,86 

3.75 

0.08 

0.04 

0.03 

0.03 

0.02 

0 

0 

0.10 

0.12 

0.07 

0.06 

0.05 

0.02 

0.02 

0.18 

0.16 

0.10 

0.09 

0.07 

0.02 

0.02 

1.83 

1.33 

1.43 

1.50 

1.40 

1.00 

1.00 

0.023 

0,060 

0.038 

0.094 

0.075 

0.026 

0.090 

0.288 

0.444 

0.518 

1.564 

1.390 

0.857 

3.747 

0.86 

1.33 

2.33 
3.42 
6.55 
9.53 
0.86 

1.151 

1.774 

2.848 

4.985 

7.940 

10.387 

4,610 

0.041 

0.167 

0.171 

0.735 

0.929 

0.419 

-0.651 

(cu  ft) 

1.811 

Crater 

0.49 

IBl 

BQ 

EQI 

0.59  1  0.59 

1.0 

0.143 

0.163 

0.375 

0.538 

0.121 

East 

■ 

0.13 

1.04 

1.45 

4.71 

2.59 

2.75 

0.07 

0.52 

0.78 

2.36 
1.29 

1.37 

0,04 

0.35 

0.48 

1.57 

0.86 

0.92 

0 

0.13 

0,04 

0.01 

0.01 

0.02 

0.13 

0.04 

0.01 

0.01 

0.02 

0 

0.13 

0.17 

0.05 

0.02 

0.03 

0.02 

0.26 

0.21 

0.06 

0.03 

0.05 

0.02 

2.0 

1.24 

1.20 

1.50 

1.67 

1.00 

0.009 

0.088 

0.039 

0.047 

0.039 

0.028 

0.088 

0.429 

0.580 

2.355 

1.441 

0.917 

1.03 

1.16 

2,20 

3.65 

8.36 

10.95 

1.117 

1.589 

2.780 

6.005 

9.801 

11.867 

0.016 

0.220 

0.170 

0.443 

0.600 

0.522 

^■1 

1.971 

H 

m 

1  0.22 

0.201 

■■ 

1 

2 

3 

1.20 

2.26 

1.97 

0.60 

1.13 

0.99 

0.03 

0.02 

0.01 

0,04 

0.03 

0.01 

1.33 

1.50 

1.00 

0.018 

0.023 

0.099 

■BSfl 

Bb 

mgmm 

1.802 

3.600 

5.3870 

0.051 

0.130 

0.838 

Tota 

L 

1.019 

IKB 

lIQByQQ 

IH 

msi 

0,308 

1  0.29 

0.665 

0.322 

Total  volume  of  crater  below  original  surface  =  U./»U 

Volume  of  hole  made  by  projectile  =  w (0.375)2  (13.5)  ,  =  5.yb4 

Volume  of  crater  above  original  surface  =  (0.07)  (tt)  (1.24)2  =  0.33« 
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Figure  60.  Description  of  Soil  Strata 
for  Test  120-57 
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ANTELOPE  UKE 

Figure  61,  Log  of  Borings  for  Tests  120-89 
and  120-97 
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of  9-inch  Blunt  Nosed  Projectile 


Figure  64c.  Plan  View  of  Soil  Nose  Taken  from  9 -inch  Blunt 
Nosed  Projectile  After  Penetration  Into  Playa 
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to  the  blunt  end.  In  Test  No.  120-57  the  paint  was  only  removed  for 
6  inches  back  from  the  blunt  end  on  the  area  that  was  on  the  bottom 
side  of  the  trajectory.  All  three  tests  were  about  5  degrees  off  the 
vertical  at  time  of  impact.  When  the  projectiles  were  removed,  soil 
was  stuck  to  the  sides,  but  the  paint  was  intact  under  the  soil.  The 
hole  remained  open  in  all  three  tests  and  none  of  the  projectiles  could 
have  been  removed  had  not  a  hole  4  feet  in  diameter  been  bored  down 
beside  the  projectile.  This  large  hole  allowed  the  projectile  and  the 
soil  nose  to  be  excavated.  The  penetration  hole  had  to  be  larger  during 
penetration  process  or  all  of  the  paint  on  the  projectile  would  have  been 
removed.  The  ultimate  angle  of  internal  friction  and  the  cohesive 
strength  of  the  soil  near  Test  Nos.  120-89  and  120-97  was  about  41 
degrees  and  1  ton  per  square  foot. 

b.  Sharp  Nosed  Tests  Into  Play  a  --  The  soil  for  these  tests  was 
similar  to  the  soil  described  for  Test  No.  120-57.  These  tests  were 
conducted  by  Sandia  Corporation  using  the  same  type  equipment  de¬ 
scribed  for  the  blunt  nosed  tests.  Soil  noses  were  not  looked  for  nor 
was  the  surface  surveyed  after  impact.  For  one  2.  2  caliber  test,  none 
of  the  paint  was  removed  immediately  behind  the  nose.  In  other 
sharp  nosed  projectile  tests,  where  the  nose  shape  was  varied,  it  is 
believed  that  the  area  of  paint  removal  immediately  behind  the  nose 
increased  as  the  nose  was  made  longer.  The  motion  time  data  for 
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three  tests  are  given  in  Figures  67,  68,  and  69.  The  depth  of  penetra¬ 
tion  listed  on  the  figures  is  a  measured  depth  where  the  penetration 
plotted  is  a  calculated  value  obtained  by  numerically  integrating  the 
deceleration  twice  with  respect  to  time.  The  deceleration  is  certainly 
not  a  linear  function  of  the  velocity  squared  at  any  instant  of  time  as  the 
Poncelet  equation  assumes.  Soil  has  been  observed  stuck  to  the  point 
in  several  of  the  tests.  The  data  presented  do  not  show  a  deceleration 
spike  upon  impact. 

By  making  the  projectile  nose  sharp,  the  surface  deceleration 
spike  must  be  highly  mitigated,  but  for  high  velocity  impact  (1000- 
1500  ft /sec),  the  afterbody  of  the  projectile  was  broken  in  tension  when 
no  deceleration  spike  was  recorded.  A  very  short  duration  deceleration 
spike  must  have  existed  in  order  to  propagate  a  compression  wave  up 
the  projectile  so  that  a  tension  wave  would  be  reflected  off  the  rear  of 
the  projectile  to  break  the  afterbody.  Generally,  the  sharp  pointed  pro¬ 
jectiles  penetrated  further  into  the  dry  lake  playa  than  blunt  nosed  pro¬ 
jectiles  when  all  other  parameters  were  held  constant. 

c.  Sharp  Nosed  Test  Into  Sand  --  Although  many  sharp  nosed 
tests  have  been  performed  by  Sandia  Corporation  into  many  different 
soils,  one  test  into  sand  is  of  particular  interest.  The  test  was  per¬ 
formed  at  Eglin  Air  Force  Base  using  a  projectile  with  a  9.  25  caliber 
point.  The  data  obtained  for  deceleration,  velocity  and  penetration 
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are  shown  in  Figure  70.  Photographs  of  the  projectile  before  and 
after  penetration  are  shown  in  Figures  71  and  72.  The  crater  and  hole 
formed  by  the  projectile  are  shown  in  Figure  73.  In  this  test,  a  soil 
nose  certainly  did  not  form  since  the  point  is  ablated.  Possibly  the 
ablation  would  not  have  happened  if  the  point  had  not  been  so  sharp. 
Knowing  that  soil  noses  do  form  on  projectiles  with  short  noses,  and 
that  long  nosed  projectiles  ablate,  it  seems  reasonable  to  assume 
there  is  some  nose  shape  where  neither  phenomena  occurs.  The  line 
of  separation  between  the  soil  and  afterbody  was  about  4  feet  back  of 
the  point  on  this  9  inch  projectile. 

No  deceleration  record  was  obtained  after  one  tenth  of  a  second 
of  time  which  corresponds  to  56  feet  of  penetration  although  the  pro¬ 
jectile  continued  in  penetration  to  a  depth  of  71  feet. 


NOSE  — 9.25  CALIBER  OGIVE 
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Figure  71.  Photograph  of  9.  25  Caliber  Projectile 
Before  Penetration  Into  Sand 


IV.  SUMMARY  AND  RESULTS  OF  EXPERIMENTAL 

OBSERVATIONS 

1.  As  shown  by  all  the  deceleration  records,  the  soil  force 
which  resists  projectile  motion  is  not  constant  and  independent  of 
velocity  as  Euler  predicted  nor  is  it  a  linear  function  of  the  velocity 
squared  term  as  the  Poncelet  and  Petry  formulas  assume. 

2.  Before  the  projectile  hits,  the  state  of  stress  in  the  half-space 
of  soil  is  as  follows: 

a  =  zy 
z 

a  =  kz7  (III- la) 

r 

T  =0 

rz 

where 

7  =  effective  unit  weight  of  soil 

and 

k  =  coefficient  of  earth  pressure  at  rest, 
k  <  1 


153 


154 


The  velocities  are  all  zero  before  impact,  thus 

V  =  V  =  V  =0  (III- lb) 

r  z 

3  8 

The  initial  state  of  stress  for  any  other  target  is  different. 

During  the  event  when  r  -♦oo,  or  when  z  oo  the  equation  for  the 
stresses  and  velocities  before  impact  must  hold.  This  would  not  be 
true  for  another  kind  of  target. 

3.  A  soil  nose  forms  on  blunt,  hemispherical  and  short  ogive 
nosed  projectiles  as  the  projectile  enters  the  soil.  This  nose  is  nearly 
conical.  On  long  ogive  nosed  projectiles,  the  metal  ablates.  The 
shape  of  nose  and  angle  of  incidence  highly  affects  the  deceleration 
spike  at  time  of  impact.  The  energy  expended  as  the  projectile 
passes  from  air  into  the  soil  seems  small  because  of  the  very  short 
duration  of  the  deceleration  spike. 

Along  the  trajectory  below  the  projectile,  there  is  no  movement 

of  the  soil.  Thus  t  is  equal  to  zero  at  the  center  line.  From  the 

rz 

Mohr  circle.  Figure  11,  it  can  be  seen  that  the  direction  of  flow  /3^ 

is  -7  -  ^  and  the  pseudo  direction  of  yield  or  flow  ^  =13  +  ^  +  (j) 

4  ^  n  s  ^ 

o  o 

R.  O.  Hultgren,  L.  J.  Thompson,  and  D.  G.  Kitzinger,  State 
of  Stress  in  an  Elastic  Earth  Target  Resting  on  a  Rigid  Half -Space. 
Sandia  Corporation,  Albuquerque,  New  Mexico,  Research  Report 
64-1734,  May  1965. 


when  is  equal  to  zero.  The  nose  on  a  blunt  nosed  projectile  as 
seen  in  the  two  and  three  dimensional  tests  has  edge  slopes  near  these 
angles.  In  Figure  74a  it  can  be  seen  that  the  equation  of  the  edge  of 
the  nose  cone  is 


(III-2a) 


for 


d  <  z  < 


R 

P 


2 


) 


and 


0  <  r  <  R 

-  -  p 


The  mass  of  the  soil  cone  does  not  change  with  time;  therefore,  M 


is 


constant. 


The  equation  for  ,  the  conical  mass  of  the  soil  carried  along 
with  a  blunt  nosed  projectile,  is 


M 

c 


R" 


tan  I 


/  rr 


(III -2b) 
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Moving 


Figure  74.  State  of  Stress  and  Velocity  at  the  Soil  Nose 
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4.  As  the  projectile  impacts,  the  soil  is  pushed  down  and 
outward.  In  the  short  interval  of  time  it  takes  the  projectile  to  pass 
through  the  interface,  the  soil  particles  go  from  zero  velocity  to  some 
high  velocity.  The  direction  of  motion  during  this  short  interval 
seems  to  follow  a  log  spiral  centered  at  the  edge  of  the  face  of  the 
blunt  nosed  projectile.  In  order  to  separate  out  the  interface  effects, 
it  looks  reasonable  to  assume  for  initial  conditions  that 

a.  The  soil  nose  has  formed, 

b.  The  direction  of  the  velocity  of  the  soil  particles  follows 
a  log  spiral  centered  on  the  edge  of  the  projectile  face. 

Below  the  log  spiral  that  has  an  initial  radius  of  the  length 
of  the  side  of  the  point,  no  motion  is  occurring.  The 
velocity  of  the  soil  particles  is  a  function  of  the  velocity 
of  the  projectile  and  position  of  the  particles  only. 

c.  The  projectile  is  moving  at  a  velocity  which  is  less 
than  the  velocity  before  impact. 

A  lower  bound  on  the  change  in  velocity  at  the  interface  can  be 
obtained  if  it  is  assumed  the  soil  offers  no  resistance  to  the  movement 
of  the  soil  nose  formed.  If  the  soil  offers  no  resistance,  momentum 
would  be  conserved  and 

M 

'^o  M  +  M  i 
P  c 
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Actually  the  target  resistance  at  time  of  impact  is  high  so  this  is  a 
poor  approximation,  and  momentum  is  not  conserved;  thus 


M 

V  < - 2 - V 

''o-  M  +  M  i 
P  c 


Or  the  velocity  just  after  impact  can  be  assumed  as 


(III -3  a) 


V 

o 


M 

P 


k  M  V. 
P  1 


(III- 3b) 


where  k  depends  on  the  nose  shape,  projectile  diameter,  target  soil, 
and  even  perhaps 

5.  After  the  projectile  has  entered  the  soil,  the  stresses  on  the 
soil  nose  can  be  found  from  the  Mohr  diagram  when  its  shape  is  known. 
From  Figure  11  and  from  Equations  11-34,  11-35,  and  11-36  it  can  be 
seen  that 


a 

z 


a  + 
s 


C  +  a 

s 


/9  V 

tan^ +  K(-j^  + 


1  +  sin  (j) ' 
cos  0 


a  = 
r 


C  +  tan  0  + 


8  V 
_ 1 

9z 


+ 


9  V 

2 

9r 


)li 


1  -  sin0 
cos  0 


T  =  0 

rz 


when  /3  = 


1 

2 
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The  shear  stress  parallel  tothe  soilnose  is  r  and  the  stress  normal  to 

n 

soil  nose  is  But  from  Equation  11-32  and  the  Mohr  circle  in 
Figure  11  it  is  seen  that 

/av  av  \ 

T  =C  +  CT  tan0  +  k(^-^  + )  (III-4) 

n  n  \  az  ar  / 

CT  =  CT 

n  s 

where  the  direction  of  s  is  /3  and  the  direction  of  n  is  /3  +  ^  +  <j), 

s  s  —  ^ 

since  both  s  and  n  are  yield  and  flow  lines. 

As  shown  in  Figure  74a  the  vertical  component  of  is 

T  =  T  cosf-j  - 

nz  n  \4  2/ 

and  the  vertical  component  of  cr  is 

n 

a  =  or  sin(—  -  —] 

nz  s  \4  2/ 

The  sum  t  +  ct  is  the  vertical  force  on  the  differential  area 
nz  nz 

element  dA  of  the  soil  nose.  Since  the  cone  is  symmetric  about  the 
centerline  z,  dA  can  be  taken  as  the  area  of  a  strip  dz  in  height  around 
the  nose  where 


or 
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dA  =2/7  CSC 


V4  2  / 


r  dr 


since 


dr 

dz 


at  the  nose 


The  total  force  on  the  nose,  F,  is  the  integral  of  the  forces  on  the 
differential  areas.  Thus 


r=0 


r  dr 


(III- 5) 


By  substitution  of  Equation  III- 4 


F  =  2;7  I 

^r=R 

Pf 

^  +  tan  (j>  cot  ~  2^ 

)  *  c  (I  - 1) 

]  rdr  (III- 6) 

•i 

'=0 

where  ct  is  as  function  of  depth  and  velocity  of  the  projectile  and 
s 

position  on  soil  nose  cone  at  any  instant  of  time. 

6.  Since  the  sides  of  the  soil  cone  seem  to  be  straight  except  at 
the  point,  the  flow  lines  must  be  parallel  to  each  other.  Then  over 
the  whole  nose  the  angle  to  the  s  line  (direction  of  flow)  is  ^  ^ 

and  the  angle  /3^  to  n  line  is  (3^  +  ^  +  <j).  If  the  flow  lines  are  parallel 
at  the  surface  of  the  cone,  the  velocity  vectors  of  the  soil  being  dis¬ 
placed  must  be  uniform  over  the  nose.  All  points  on  the  nose  move 
down  the  same  distance  during  the  short  time  interval  At. 
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The  distance  the  soil  particles  naove  is  U  in  the  direction 
^  ~  2'  ^  seen  in  Figure  74b,  c,  and  d, 

that 

V 

V  = - 7 - 7-  in  the  direction  B  -  45  -  ^  (111-7) 

n  9\  2 

2  oos(j-2) 


(111-8) 


As  the  blunt  nose  projectile  travels  through  soil,  there  seems  to 
be  little  friction  on  the  sides  of  the  projectile.  An  open  cavity  sur¬ 
rounds  the  projectile  except  at  the  nose.  Along  the  edge  of  the  cavity 


T  =0 
ns 

(111-9) 

a  =  0 
n 


and 


V  =  0 
n 


where  the  direction  s  is  along  the  cavity  and  the  direction  n  is  normal 


to  the  surface  of  the  cavity. 
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A  reasonable  assumption  for  the  equation  of  the  cavity  seems 

to  be 

r  =  R  +  (d  -  z)  k  (III- 10) 

P 

where  z  <  d,  r  >  R^  and  k  is  constant  and  between  and 

zero. 

To  obtain  stability  of  the  projectile,  its  length  must  be  several 
times  its  diameter  so  that  it  does  not  tumble  in  the  cavity.  Possibly 
the  length  to  diameter  ratio  must  be  at  least  7  or  8  to  1  to  obtain  the 
desired  stability  for  large  projectiles. 

8.  For  the  penetration  hole  to  remain  open,  the  soil  must  arch 
around  the  cavity  created  by  the  projectile.  This  only  happens  if  the 
soil  has  some  cohesive  strength.  After  the  projectile  has  stopped, 
the  hole  remaining  open  in  the  soil  is  somewhat  smaller  than  the 
projectile.  Therefore,  on  the  surface  of  the  projectile,  there  are  normal 
forces  which  cause  frictional  forces  that  resist  removal  of  the 
projectile;  but  these  forces  do  not  seem  to  exist  during  penetration, 
since,  in  the  case  of  the  blunt  nosed  projectile,  the  paint  is  not 
removed  even  though  soil  is  stuck  to  the  projectile. 

9.  The  stress  a  and  t  are  zero  at  the  surface.  Thus  a 

z  rz  z 

is  the  minimum  principal  stress  and  equal  to  cr^.  The  cracks  at  the 
surface  show  that  the  radial  normal  stress  has  exceeded  the  strength 
of  the  soil  in  tension  or  that  ct  is  equal  to  zero. 


10.  The  volume  of  the  soil  raised  at  the  surface  is  about  the 
same  as  the  volume  of  the  hole  created  by  the  projectile. 

The  heave  of  soil  at  the  surface  shows  that  the  assumption  that 
the  density  is  constant  is  not  unrealistic. 

11.  The  motion  of  the  target  in  the  two-dimensional  tests,  as 
weU  as  the  length  of  surface  cracks  and  surface  deformations  in  the 
full  scale  tests  indicate  that  the  phenomenon  of  penetration  is  primarily 
one  of  shear  deformation,  if  the  soil  density  is  above  critical  density. 
There  seems  to  be  a  shear  front  that  travels  with  the  projectile.  The 
equation  of  the  leading  edge  of  the  front  seems  to  be  a  log  spiral  in 

the  case  of  a  blunt  nosed  projectile  penetrating  a  half  space.  For 
other  types  of  targets,  the  equation  of  the  shear  front  probably  depends 
on  the  boundary  constraints. 

The  motion  of  the  particles  in  the  two-dimensional  tests  indi¬ 
cates  that  the  shear  front  translates  with  the  projectile  and  that  the 
velocity  of  the  target  particles  is  a  function  of  the  position  and 
velocity  of  the  projectile. 

For  pure  translation  of  the  velocity  field,  which  is  a  function  of 
the  velocity  of  the  projectile,  the  components  of  the  velocity  field  can 


.be  written 
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where  V  can  be  treated  as  a  constant  over  the  whole  field  at  an  instant 
P 

of  time. 

The  acceleration  terms  for  the  field  particles  are  then 


d  V 

a  V 

a  V 

dz 

r 

r 

r 

dt 

dr 

dt  ^ 

az 

dt 

and 


dV 

a  V 

^  + 

a  V 

dz 

z 

z 

z 

dt 

ar 

dt 

az 

dt 

dV  a  V 

- ^  =  V  - ^ 

dt  p  az 


(III- 12) 


12.  It  is  observed  that  there  is  a  discontinuity  in  the  displacement 
and  velocity  fields  below  the  point  and  to  the  side  of  the  soil  that  is 
moved.  This  discontinuity  defines  the  shear  front.  On  the  edge  of  this 
discontinuity,  the  velocity  components  are  zero. 
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13.  It  can  be  seen  from  the  two-dimensional  test  photographs 
that  at  any  instant  of  time  the  displacement  of  the  particles  is  very 
similar  to  the  displacements  initially  approximated  by  Prandtl  for  the 
indentation  of  an  infinitely  long  footing  or  punch  into  a  half  space  of 
rigid-plastic  weightless  material.  There  is  a  zone  of  material  below 
the  projectile  that  has  not  been  influenced  by  the  presence  of  the  pro¬ 
jectile.  If  the  movement  of  the  target  particles  is  along  a  log  spiral 
centered  on  the  edge  of  the  projectile,  then  there  can  be  no  velocity 
component  perpendicular  to  the  log  spiral. 

K  the  soil  does  move  at  each  instant  of  time  along  a  log  spiral 
centered  on  the  edge  of  the  projectile,  the  velocity  field  can  be  calcula 
ted  as  follows. 


Let  X  =  A  q  e^  ^  where  x  is  a  new  variable, 


A  is  a  constant  =  e 


(i-i) 


tan  0 


and  u  is  a  new  variable  such  that 


w  +  T  =  j3  +  0- 


From  Figure  75  it  can  be  seen  that 


r  =  R  +  X  sinu  =R  +sinu  (Aqe 
P  P 


u  tan  0 


u  tan  0 


) 


(III- 13) 


z  =  d  +  X  cosoj  =  d  +  cosu  (Aqe 


(III- 14) 


Direction  of  Flow  Lines  at  an  Instant  of  Time 


By  taking  the  partial  derivative  of  Equations  111-13  and  111-14  with 
respect  to  r,  the  two  following  equations  are  obtained. 
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8r 


8u) 


V  OCJ 

+  (tanp  +  cotu)  q 


A  sinw  e 


u)  tan  0 


(III- 15) 


and 


8r 


+  (tan^  -  tanuj)  ‘I  ^  ^ 


(III- 16) 


Solving  simultaneously  Equations  III- 15  and  III- 16 


8cj 

8r 


cos  CJ 


q  Ae 


tan  <j) 


8r 


cosu)  tan0  ^  sinco  I 
Ae^tan^  ^^cotan(^J 


(III- 17) 


Taking  the  partial  derivative  of  Equations  III- 13  and  III- 14  with  respect 


to  z 


|a+(tan^  +  cotu,)q^ 


=  0 


and 


8q  ^  ,  ,  V  8u 

+  (tanp  -  tanu)  q 
8z  8z 


_ 1 _ 

.  tan0 

A  cosu  e 
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Using  elimination 


du 

9z 


sincj 


q  A  e 


uj  tan  (j) 


9q  _  sincj  tan  ^  ^  cos  u 

~  .  Lj  tan  (j)  .  u  tan  (f) 

A  e  A  e 


Since 


n 


13  =  (J  +  -  </> 


98  9to  ,  9^  9u 

=  TT-  and  =  -r— 

9r  9r  9z  9z 


The  continuity  equation  with  p  =  0  is 


9  V  V  9  V 
r +_r  ^  2 


9r 


9z 


=  0 


where  V  and  V  are  functions  of  q  and  co  at  any  instant  of  time, 

X*  2i 


9V,  9V^  V  av„  9V 

_ r  ^  _ r  ^  ^ _ z  ^  _ z  ^ 

9q  9r  9  to  9r  r  9q  9z  9  to  9z 


=  0 


But 


V 

r 


V 

z 


=  V  sinjS  =  V  sin^to  +  ^  -  ^  =  ^  cos|(o  - 


=  V  cos/S  =  V  COSTCO  +  ^  "V  sin^to  -  ^ 


(El- 18) 


(11-19) 

thus 


(IE- 19) 
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thus 


(III-20) 


By  substituting  Equations  III- 17,  III- 18  and  III- 20  into  the  continuity 
Equation,  11-19,  the  following  equation  is  obtained. 


cos(w  -  +  [  cos(w  -  <^)|^  -  V  sin(w  -  (/.)j  ^  cos(w  -  0) 

-  sinu  -  0  -  {  sin(u  -  </>)  |^  +  V  cos(cj  -  0)]  =  0 


or 


av  cos(u)  -  sin  CO  -  cos  cj  tan  0)  -  sin(u  -  ^Xcosoj  +  sinu  tan^ 


+ 


9V  f  cos(  u  -  <^)(cosco)+  sin(u)  -  0)(sinuj) 


9tj 


q  A  e 


utan  <j) 


+  V 


rcos(uj  -  0)  sin(cL)  -  (^)cosuj  ^  cos('j 

I  r  »  wtan0 

^  q  Ae 


170 


which /reduces  to 


9V  cos  0  ^,/cos((j  -  <^) , 

9u  .  to  tan  0  \  r 

q  A  e 


sin(/)  V  Q 

.  k)tan(^  ) 
q  Ae  ' 


or 


1  9V  ,  .  ,  ,  COS(u3  - 

—  -7^  +  tan  q  + - -jr 

V  9uj  COS  9 


.  wtan^ 

qAe _ ^  Q 


R  +  Aqe 
P 


utan  (j) 


(III-21) 


Since 


and 


^  ,  .  ijtan^  . 

r  =  R  +  Aqe  sinw 

P 


-^=  Aqe^*^''*^cosu;  + Aqe"^^^”*^  sincotan.^. 
9w 


then 


1  9v  ^  ^  ,  1  ar  _  o 

V  du)  ^  r  au 


(III-22) 


Integrating  with  respect  to  u,  it  is  found  that 
in  V  +  utan(^»  +  i  n  r  +  f(q)  =  0  ,  or 


V  =  e  ^  e  which  reduces  to 


V 


-utan(/>  -f(q) 
e  e 


(III-23) 
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At  the 

where 

Then 


Then 


soil  nose 


for 


0  <  q  < 


V  = 


V  h 
P _ 

^  utanc^ 
Be  r 


(III-24) 


(III-25) 
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where  h  =  R  -  q  sin^-|  -  ^ 


B  =  2  A  cos 


(f  - 1) 


■O  A  Cjtan^  . 

r  =  R  +  Aqe  sinu 

P 


V 


V  h  cos  (u  - 


r  -n  w  ^ 
Be  r 


V 

z 


V  h  sincj  - 

_£ _ 

T_,  u  tan  0 
Be 


<!> 


r 


(III- 26) 


(m-27) 


when  the  velocity  of  the  projectile  V  is  equal  to  the  velocity  after 

P 

impact  V 


rr  ,  A  .  ^  n 

4  2  —  —  2 


By  substituting  '\^  for  into  Equation  III-21  and  allowing  co  to 
vary  between  -nlA+  0/2  and  /t/2,  the  initial  velocity  field  can  be  estab¬ 
lished.  When  Vp  =  0  the  projectile  has  stopped.  The  assumed  log  spiral 
deformation  and  velocity  field  require  that  and  be  zero  at  this 
time. 

It  is  also  observed  from  the  two-dimensional  studies  that  when  to 
is  between  rr/2  +  0  and  njA  +  30/2  that  the  soil  leaves  the  log  spiral 
motion  and  moves  upward  and  outward  in  an  almost  straight  line  parallel 
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to  the  side  of  the  cavity.  The  above  limits  for  o)  correspond  to  jS  by 
being  between  n  and  SnlA  +  0/2. 

If  the  soil  does  move  along  lines  in  which  jS  is  constant,  then  the 
velocity  field  above  the  log  spiral  area  can  be  calculated  as  follows. 

Let 

r  =  R  +  Q  cos  (0  -  O')  +  W  sin  a  (III-28) 

z  =  d  -  Q  sin  (0  -  of)  -  W  cos  a  (111-29) 

As  shown  in  Figure  76,  V  and  V  then  would  be  functions  of  Q  and  W 
and  at  an  instant  of  time.  Thus 

^  ^  9V  3Q  ^  9V  8W 
9r  9Q  9r  9W  9r 

and 

9z  9Q  9z  9W  9z 

Also 

V  =  V  sin  S  =  V  sin  {n  -  a)  =  V  sin  a 
r 

V  =  V  cos  jS  =  V  cos  (n  -  a)  =  -V  cos  a 
z 


where  a  is  a  constant  between  0  and  niA  -  0/2. 
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The  continuity  equation  (11-19)  can  be  rewritten 


sin  O' 


/dV  dQ  dV  8W\  V  sin 
VaQ  ar  aw  ar  /  r 


a 


cos 


I 


/av  aQ  ^ 
"  aQ  az 


^  9W\  ^  Q 

aw  az  j 


or 


aQ\ar 


cot 


9Q\  ,  9v/aw 
"  az/  aw\Br 


-  cot 


V 

az ) 


r 


(III- 30) 


The  partial  derivative,  dQ/dr,  aQ/8z,  dW/dr,  aw/az  can  be  evaluated 
as  follows.  Taking  the  partial  derivative  of  Equations  III-28  and  III-29 
with  respect  to  r. 


ar 


1  =  cos  {<b  -  a) 


a_Q 

ar 


+  sin  a 


dr 


az  -  .  .aQ  aw 

—  =  0  =  -sin  (d  -  a)~  -  cos  a  — 
dr  dr  dr 


and  solving  simultaneously 


^  _ 

cos  a 

ar 

cos 

aw  _ 

sin  -  a) 

a  r 

cos 

(III- 31) 
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Taking  the  partial  derivative  of  Equations  III-28  and  III-29  with  respect 
to  z. 


Br 
9  z 


0  =  cos  ^ 


m 

Bz 


9  z 
Bz 


1  =  -sin  (<^  - 


vBQ 

Q-lr—  -  COS  a 
9  Z 


m 

B  z 


and  solving  simultaneously 


9Q  _  sin  a 
Bz  cos  0 

9W  _  cos  -  a 
Bz  cos  <f) 


(111-32) 


Substitution  of  Equations  III-31  and  111-32  into  Equation  111-30  yields 

9V/cos_a  cos  a  sin 

9Q\cos  (j)  sin  a  cos  (j)/ 

9V  /sin((^  -  a)  _  cos  or  cos((/)  -  g)\  ^  V  _  ^ 

9W\  cos  4  sin  a  cos  0  /  r 


which  reduces  to 


1  9  V 

sin  a  9W 


or 


1  ^  sin  g 
V  m 


r 


0 


(111-33) 
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Since 


r  =  Rp  +  Q  cos{(f)  -  0')+  W  sin  a 


and 


9r 

Fw 


=  sin  a 


then  Equation  111-33  can  be  rewritten 


1^1^  =  0 
V  aw  r  8W 


(III- 34) 


or 


Integrating  Equation  III-34  with  respect  to  W  it  is  found  that 


/n  V  +  /n  r  +  f(Q)  =  0 


-in  r  -f(Q)  e 
V  =  e  *6  =  - 


-f(Q) 


(III -35) 


when  (o  =  n/2  +  ^  -  a.  Equation  III-25  can  be  written 


+  ^  -  c)  tan  0 


V  = 


Vp  [Rp  -  q  sin(f  -  |)] 


2rA(cos|-|) 


(III- 36) 


where 
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Equation  III-36  allows  the  evaluation  of  e  .  Wheno)  -  ir/2  +  cf  -  u 
Equation  III- 36  must  equal  Equation  III-  35,  thus 


-f{Q)  e 


_(|  +  ^  tan  0  .  (rr  (i)W 

'  V  Ir  -  qsin(^-j^ 


2rA  cos 


(f-l) 


-f(Q) 


-(1+ ^  -  <i)tan 

=  vlR-q  810(5-2)) 


2A  cos 


(f-l) 


(111-37) 


Substituting  Equation  III- 37  into  Equation  III- 3 5,  the  velocity  is 
found  for  Area  II  where  the  soil  is  moving  parallel  to  the  cavity 


■(■|+  (^  -  a)  tan  ^ 


^11  = 


^  ^  .  In  (l>\] 

V^[Rp-qsin(--3)| 

2rA  cos{^^  - 


(III- 38) 


CHAPTER  IV 


SUMMARY  OF  GOVERNING 
DIFFERENTIAL  EQUATIONS  AND 
REQUIRED  BOUNDARY  CONDITIONS 

I.  SUMMATION  OF  EQUATIONS 

Governing  equations  from  continuum  mechanics  are  as  follows: 

dV 

1.  (M  +  M  )  g  -  F  =  (M  +  M  )  (Eq.  I-l) 

p  c  °  p  c  dt 

or 


dV 

(M  +M)g-F  =  (M  +M)V 
p  c  °  p  c  p  dz 


which  can  be  integrated  from  a  small  increment  of  time  after  impact 
when  the  velocity  of  the  projectile  is  equal  to  V^  and  the  depth  z  is 
equal  to  zero  to  the  time  when  the  velocity  of  the  projectile  is  zero 


and  the  depth  of  penetration  is  a  maximum. 
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2  F  =  I  P  da 


(Eq,  11-20) 


Oa  9t  /  ct  -  cr  \  /dV 
r  rz  /  r  v  \  / 


9r  9z 


9  V  9  V  \ 

^1+V  -^+V  -^)  p  (Eq.  II-IO) 
9t  r  9r  z  9z  / 


9t  9cr  T 
^  rz  z  rz  , 


/9V  9V 

[ - -  +  V  - -  + 

\  9t  r  9r 


9V 

(Eq-n-ii) 


9V  V  9  V 

- -  +  —  + 

9r  r  9z 


0  (Eq.  11-18) 


Constitutive  equations  are  as  follows: 


6.  CT  =  CT  +IC+CT  tan (^  +  K I  — —  +  — ^ 
z  s  I  s  '  9z  9r 


/ 9  ^r  ^  ^  ^z  \  2  cos/3  sin(j3  +  ^) 

V  9z  9r  / 1  I  cos  ^ 


(Eq.  11-34) 


^  ^  \  2  sing  cos(3  +  ^)1 

9r  /  J  I  cos  6  1 


-  lC+^stan0  +  K(^— +  -^ 


cos  (j) 


(Eq.  11-35) 


(  ^  ^  ^  [  cos (2  s  +  ^1 

\  9z  9r  I  cos  0  ] 


(Eq.  11-36) 


9.  tan^  =  ~  (Eq.  11-39) 
z 


10. 


p=  constant  (Eq.  11-42) 
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(Eq.  in-6) 


II.  BOUNDARY  AND  INITIAL  CONDITIONS 
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There  are  13  equations  with  13  unknowns.  The  initial  velocity- 
condition  when  d  =  0,  and  for  the  soil  outside  the  nose  and 

above  the  shear  front  discontinuity  is 


V  = 


V  h 
P 


_  tjtan  A 
Be  ^ 


(Eq.  in-21) 


and 


V  =  V  sin 
r 

V  =  V  cos  B 
z 

from  Equation  n- 39.  For  the  soil  inside  the  soil  nose,  and 

V  =  0.  The  soil  nose  is  described  by  the  equation 
r 

R  -  r 

Z  =  - 2 -  +  d  {Eq.  III-2a) 

tan  ^45°  -  ^ 


The  final  velocity  condition  is  when  d  =  d  and  V  =  0. 

XXXa-X.  ^ 


The  state  of  stress  and  velocity  before  impact  is; 


CT  =  z 
z 


a  =  z 
r 

r  =0 
rz 


V  =  V  =  V  =  0 
r  z 


V  (Eq.  III-D 


F or  all  time  during  the  event  the  botindary  conditions  for 
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a.  the  surface  z  =  0 


and 

V  ,  V  and  cr  are  unknown, 
r  z  r 

b.  Along  the  cavity  described  by 


r  =  R  +  (d  -  z)  k  (Eq.  Ill- 10) 

\  ‘  \  -  °  “1-9) 

and 

V  and  <j  are  unknown  on  this  boundary. 

So 


c. 

velocity 


Along  the  nose  when  the  projectile  is  at  depth  d  and 
Vp,  the  velocity  is  uniform  and  given  by 

V  h 


V  = 


P 

(j  tan  6 
Be  ^ 


(Eq.  Ill- 21) 


where 


4  2 


0)  = 


and 


184 


=  vcos(a-|) 


is  luiknown  on  this  boxmdary. 

d.  There  is  a  shear  front  which  travels  with  the  projectile;  the 
equation  for  the  shear  front  is  unknown  theoretically,  but  outside  this 
front  where  =  0,  the  boundary  stresses  along  the  front  can  be 

obtained  from  Equations  11-34,  11-35,  11-36,  and  11-41. 


a 

z 


=  cr  +  (C  +  (j  tan  0 
s  s 


2  cos  j8  sin  ()6  +  <j)) 
cos  ^ 


(C  +  CT 

s 


tan  4>) 


2  sin  /3  cos  (/3  +  <!>) 
cos  (j> 


o  ,  =  iC  +  a 

yj  S  S 


tan 


sin  (j)  -  1 
cos  0 


-(C  +  tan  (j)) 


cos  {2(3+  <l>) 
cos  <f> 


V  =  V  =  0 
r  z 

The  above  13  equations,  the  one  initial  condition,  and  the  four 
boimdary  conditions  completely  describe  the  problem. 

Outside  the  boundaries  where  the  soil  is  at  rest,  i.  e. ,  when  r  or 
z  is  large,  the  initial  conditions  before  impact  prevail. 
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m.  COMBINATION  OF  EQUATIONS 


By  differentiating  Equation  11-35  with  respect  to  r,  differentiating 
Equation  11-34  with  respect  to  z,  and  differentiating  Equation  11-36 
once  with  respect  to  r  and  once  with  respect  to  z,  the  following  equa¬ 
tions  are  obtained. 


^*^r  ^"^s  2  sin  ^  cos  0  +  <^) 


V  a^  V 

r  .  z 


ar  dr 


cos  ^ 


Ida 

tan(/)  +  Kl  + 


araz  ■  ^2 
ar 


2  cos  (2j3  +  <^) 
cos  ([> 

(Eq.  IV- 1) 


/  a  V  a  V 

C  +  tan^  +  k(  ^ 


s 


az  ar 


)i© 


^^z  _  ^°^s  ^  2  cos  ^  sinQ  +^) 
az  ~  az  cos  (f) 


da  /a^v  a^v 

'  az 


+ 


2  cos 


(20  +  0)  f 

s  I 


cos 


/  a  a 

C  +  <T^tan.>+K(^-3^  +  -^ 


(Eq.  IV- 2) 


dr 


rz  cos  (2^  +  (f>) 


dz 


cos  <f) 


aa  fa\  a\ 

tan  ^  ^  +  k1 - ^  + 


az 


a  a 

r  z> 


2  sin  (2^  +(^) 
cos  0 


'av  av 


C  +  CT  tan  0  +  K  .  - 
s  ^  \  a 


+ 


a 

r  /  J  z 


(Eq.  IV- 3) 
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9T 


rz  _  cos(2^+^ 


9r 


cos  <f> 


tan  (j) 


da 
_ s 

9r 


+  K 


'  9^  V 
_ r 

9r  9z 


S')] 


+  2  sin (2^+^  ta.n<j>  +  K{ 

cos  0  I  S  ' 


9  V  9  V 


9z 


+ 


^)1 

9r  /J 


93 

9r 


(Eq.  IV- 4) 


Subtracting  Equation  11-40  from  Equation  n-34  and  dividing  by 
(-r)  the  following  equation  is  obtained. 


\  r  /  r  I  £ 


9z  9r 


)]l^ 


cos  ^ 


(Eq.  rV-5) 


Substituting  Equation  IV- 1,  Equation  IV- 3  and  Equation  rV-5  into 
Equation  II- 10,  the  following  equation  of  motion  is  obtained. 


9ct 


9(7 


/9^  V  9^  V 

-^  +  tan  ^  ^  +  k( 

9r  cos  (j)  I  9r  \  9r  9^  ^  j 


*  ^  +  ‘'s  +  -3^)]  H 


„  r  3ff  /8^V  8^V 

+  tan*^  +k'  “■ 

’  dz 


+ 


COS 


9z 


9r  9z 


)1 


2  sin(2^  + 
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1  -  sin  (2^  + 
r  cos  <j) 
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(Eq.  IV- 6) 


Substituting  Equation  IV- 4,  Equation  IV- 2,  Equation  11-36  into 
Equation  II- 11,  the  other  equation  of  motion  is  obtained. 
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By  substituting  Equations  III- 2b  and  111-6  into  Equation  I-l,  the 
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Eq.  IV- 8) 


r  dr 


substituting  into  Equation  11-18  to  obtain  Equation  11-42 


av  V  av 

— £  +  — £  4-  — - 

ar  r  a  z 


0 


(Eq.  11-19) 


Taking  the  five  equations/  11-38,  11-19',  IV-8,  IV-6,  and  IV- 7,  it 

can  be  seen  that  the  five  unknowns  are  a  ,  V  ,  V  ,  B  and  V  :  and  the 

s  r  z  p 

three  independent  variables  are  r,  z  and  t.  B  can  be  eliminated,  leaving 
four  unknowns :  This  leaves  the  equation  of  motion  of  the  projectile 
which  requires  the  initial  condition  for  a  solution,  the  two  second 
order  equations  of  motion  which  require  two  boundary  conditions,  and 
the  continuity  condition  which  requires  one  boundary  condition. 

Theoretically  no  assumptions  must  be  made  concerning  the 
velocity  field  in  order  to  obtain  a  solution,  but  a  closed  form  solution 
to  four  simultaneous  partial  differential  equations,  two  of  which  are 
second  order  and  two  of  which  are  first  order,  is  difficult  to  obtain. 
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It  may  be  possible  to  obtain  a  numerical  solution  to  the  governing 

3  8 

differential  equations  using  the  “marker  and  cell"  technique  on  a 
high  speed  computer.  This  technique  has  been  successfully  used  for 
plane  strain  incompressible  viscous  fluid  problems,  and  has  just 
recently  been  extended  to  axially  symmetric  problems. 


Francis  H.  Harlow  and  J.  Eddie  Welch,  "Numerical  Study  of 
Large  Amplitude  Free  Surface  Motions."  This  paper  has  been  sub¬ 
mitted  for  publication  in  The  Physics  of  Fluids,  authors  are  employed 
at  Los  Alamos  Scientific  Laboratory,  Los  Alamos,  New  Mexico, 


CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDATIONS 

Analysis  and  experimental  observations  support  the  following 
conclusions  regarding  axially  symmetric  vertical  penetration  into  a 
half  space  of  particulate  medium: 

1.  The  Euler  and  Poncelet  assumptions  concerning  the  soil 
force  resisting  motion  do  not  seem  to  be  valid. 

2.  Equations  can  be  written  to  describe  the  whole  process  of 
penetration.  The  experiments  seem  to  validate  those 
equations.  The  boundary  and  initial  conditions  for  the 
problem  can  be  delineated. 

3.  The  process  of  penetration  is  primarily  one  of  shear 
deformation  of  the  soil  when  it  is  in  a  dense  condition. 

4.  Approximations  can  be  made  for  the  velocity  and  deforma¬ 
tion  fields  that  will  contribute  to  the  solution  of  the  various 
vertical  penetration  problems. 

It  is  recommended  that  these  studies  should  be  extended  in  the 
following  ways : 

1.  A  numerical  solution  to  the  governing  differential  equations 
should  be  attempted  using  the  "marker  and  cell"  technique 
on  a  high  speed  computer. 
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2.  The  photo -plastic  two  dimensional  studies  should  be  extended, 
using  more  sophisticated  equipment  to  evaluate  the  effects 

of  boundary  conditions. 

3.  Three  dimensional  laboratory  studies  should  be  continued 
using  the  instrumentation  developed  by  Division  9226  of 
Sandia  Corporation  in  order  to  study  vertical  penetration 
into  real  materials. 

4.  The  constant  velocity  tests  into  steel  needle  bearing  rollers 
were  most  helpful  in  reaching  an  understanding  of  the  pen¬ 
etration  process.  This  work  should  be  continued  using  higher 
velocities,  various  penetration  attitudes,  different  boundary 
conditions  and  possibly  even  with  the  use  of  a  semi -flexible 
penetrator.  The  effect  of  cohesion  could  also  be  studied  by 
filling  the  voids  with  a  viscous  cohesive  material. 
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